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Abstract. Fast Ignition Inertial Confinement Fusion is a variant of inertial fusion 
in which DT fuel is first compressed to high density and then ignited by a relativistic 
electron beam generated by a fast (< 20 ps) ultra-intense laser pulse, which is usually 
brought in to the dense plasma via the inclusion of a re-entrant cone. The transport 
*yv . of this beam from the cone apex into the dense fuel is a critical part of this scheme, as 

it can strongly influence the overall energetics. Here we review progress in the theory 
and numerical simulation of fast electron transport in the context of Fast Ignition. 
Important aspects of the basic plasma physics, descriptions of the numerical methods 
used, a review of ignition-scale simulations, and a survey of schemes for controlling 
the propagation of fast electrons are included. Considerable progress has taken place 
in this area, but the development of a robust, high-gain FI 'point design' is still an 
ongoing challenge. 



* 



1. Introduction 

Since its proposal by Tabak and co-workers (Tabak et al. 1994) in 1994 the concept 
of Fast Ignition (FI) Inertial Confinement Fusion (ICF) has attracted considerable 
attention (Tabak et al. 2005). This advanced ICF concept is appealing because of 
its ability to achieve high energy gains (G > 100) while reducing both the total laser 
energy and the hydrodynamic demands on the fuel assembly. One of the new challenges 
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in this concept is the need to efficiently couple the ignitor pulse energy via the relativistic 
(fast) electrons to a hot spot in the compressed fuel. Within this, there are two parts - 
the absorption of laser light into fast electrons, and then the propagation and stopping 
of the fast electrons. This review is concerned with the latter of these, i.e. fast electron 
transport (FET). 

The fast electron transport aspect of FI is challenging for at least three reasons. 
Firstly there is an issue that would exist even if fast electron propagation were purely 
ballistic. The size of the hot spot is comparable to the size of the fast electron source (i.e. 
the laser spot), but the two are separated by a distance which is several times their size. 
Therefore any appreciable angular spread in the fast electrons must either be mitigated 
or controlled, as a reduction in the coupling efficiency will otherwise occur. Secondly 
there is the possibility that various instabilities might disrupt the beam propagation 
which in turn would impair the coupling efficiency. Thirdly, any solution to the first 
and second problem must be compatible with achievable fuel assemblies and the fast 
electron parameters required to achieve stopping in the hot spot. Yet another problem is 
the source characteristics as a function of the laser parameters: currently it would appear 
that the fast electron energy spectrum is too hard to allow for all the fast electrons to 
be deposited in an ideal hot spot. 

Ultimately it is hoped that an overarching solution to these problems can be found 
which is still attractive and feasible, i.e. a 'point design' for Fast Ignition. Currently it is 
not possible to build an ignition-scale facility purely for the purposes of investigating the 
feasbility of FI or solving the problems associated with fast electron transport by purely 
iterative empirical methods. Therefore detailed numerical simulation and a thorough 
understanding of the underlying theory are essential parts of realizing FI. Hence the 
importance of the subject matter of this review. 

In this review of the area, we will cover the following aspects of the theory of fast 
electron transport in FI: 

(i) Basic Physics: The fundamental physical processes including scattering and 
stopping of fast electrons, the role of resistively generated fields, and key beam- 
plasma instabilities and phenomena. 

(ii) Simulation Methods: The different simulation methods that have been applied to 
this problem and their relative strengths and weaknesses. 

(iii) Review of Ignition- Scale Calculations: A review of simulation studies of the full- 
scale problem, and how this has informed the overall view of the current challenges 
that FI is facing. 

(iv) Concepts for Controlling Transport: A survey of the various ideas that have been 
proposed to overcome the limits on the coupling efficiency that are imposed by 
realistic fast electron divergence angles, namely: fast electron self-collimation by 
resistively-generated magnetic fields (due to beam profile or resistivity gradients), 
electrostatic confinement by a vacuum gap (double-cone target), and imposed axial 
magnetic fields. 
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(v) Prospects for a Point Design: What future FET studies will have to address in 
order to move closer to a FI point design. 

In addition, we provide a very brief precis of the requirements for fast electron heating 
to reach ignition in section 2. 

Our review will draw attention to the considerable effort that has gone into both 
understanding the fundamental aspects of the problem and developing numerical tools 
that are suitable for studying fast electron transport. The calculations that have been 
performed under conditions close to full-scale FI so far clearly show that the scheme 
must be adapted in some way given realistic fast electron beam parameters. Our review 
also indicates there are potentially feasible ways to 'control' fast electron transport and 
thus achieve a viable point design. 

2. Ignition via Rapid Heating of Compressed Fuel 

So that the fast electron transport problem can be put in context, we briefly summarize 
the objectives that must be achieved in order to obtain ignition and gain from the 
rapid heating of a particular region of highly compressed DT fuel. Fast ignition is 
an isochoric mode of ignition, where a region of fuel of relatively constant density is 
heated to much higher temperatures and pressures. It contrasts with isobaric ignition 
modes, such as central hot-spot ignition. The requirements for FI are determined by 
hydrodynamics and burn physics. Estimates of the optimal parameters that minimize 
the ignitor pulse energy have been obtained by Atzeni using both analytic calculation 
and 2D hydrodynamic simulations (Atzeni 1999). The resulting optimal fast electron 
energy (E ign ), fast electron intensity (Ii gn ), pulse duration (t ign ), and hot spot radius 
(r hs )&re: 

E ign = 140pr o 85 kJ, (1) 

I ign = 2.4 x 10 19 p?oo 5 W cm" 2 , (2) 

t ign = 54pr °o 85 PS, (3) 

r hs = 60pr °o 9 Vm. (4) 

where pioo = p/100 g cm -3 . A typical FI scenario will involve the assembly of a quasi- 
spherical DT mass reaching peak densities in the range 300 < p < 1000 g cm -3 . 
Assuming that a re-entrant cone-guided FI scheme is being followed, the distance 
between the tip of the cone and the geometric centre of the DT mass (the 'stand- 
off' distance) is typically 100 pm. In FI schemes that employ 'hole-boring' to create a 
path for the ignitor pulse, there will still be a substantial stand-off distance of at least 
100 pm. The DT density around the cone tip depends on the detailed hydrodynamics of 
the fuel assembly scheme, but is usually on the order of at least a few g cm -3 . The fuel 
temperature at stagnation is usually around 200-300 eV. The ignition requirements were 
generalized in Ref. (Atzeni et al. 2007) to include effects like r^s exceeding the optimal 
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value and fast electrons not fully stopping in the DT fuel. More recently, ignition 
requirements based on realistic, PIC-based fast electron sources have been found (Bellei 
et al. 2012). 

The objective of fast electron transport theory is to ensure that a hot spot is 
produced within the constraints of equations (1-4) given conditions that do not differ 
too greatly from those outlined in the preceding paragraph. As Direct Drive ICF with 
central hot spot ignition may be possible with total laser energies of 1-2 MJ, and 
advanced Indirect Drive ICF may be feasible with similar total laser energy, the FI 
concepts ideally aims to operate using not much more than 100 kJ of ignitor pulse 
energy. 

3. Basic Physics of Fast Electron Transport in FI 

We now look at the basic physical phenomena that affect the propagation of the fast 
electron beam from the source to the compressed core. We will assume that the 
characteristics of the fast electron beam (FEB) at the source are given, and concentrate 
on the theoretical models describing propagation of fast electrons. In this section, 
phenomena are considered in isolation, concentrating on the fundamental equations 
and models of each phenomenon. Naturally, the interaction of these phenomena does 
occur, is rather complex, and requires use of simulation codes for quantitative prediction. 
However these fundamental elements are the 'building blocks' of fast electron transport 
theory, and are essential for understanding the simulation codes. 

3.1. Fast Electron Parameters 

The physics of the absorption of high power, high intensity laser light into the plasma 
and generation of the fast electron beam is a whole topic in its own right and is covered 
in detail elsewhere in this special issue. Such details will not be covered, and we limit 
ourselves to a few general remarks. Broadly speaking, the fast electrons are injected 
from the laser-plasma interaction region (e.g. inner cone surface) towards the compressed 
fuel core with a broad distribution of energies and a significant degree of anisotropy. A 
simple model that is sometimes used in transport calculations is, 

f inj (E, 6) oc eM-E/E) exp ( -Jp- J • (5) 

E = m e c 2 [^i — 1) is the fast electron kinetic energy, and 9 the angle between v and the 
nominal direction of beam propagation. 

The mean fast electron energy is often taken to be close to the ponderomotive 
potential energy, U pond = [1 + 0.73I 18 (X L )IJ 1/2 _ 1 (/ 18 = J/1 18 W cm" 2 ), in such 
simple models of the fast electron distribution. Fast electron transport calculations 
that aspire to have good predictive capability need to include either a self-consistent, 
laser-generated fast electron source or take a detailed fast electron distribution from a 
PIC LPI calculation. The ponderomotive scaling does provide a rough indication of 
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the expected fast electron mean energy or temperature with intensity and wavelength, 
but the detailed scalings are still very much an open topic of research. Detailed PIC 
simulations often show an energy distribution that is considerably more complex than 
the exponential in Eq. (5). 

The angular spread of the fast electrons has no simple or clean characterization 
either, all the more so since it is expected to depend on the electron energy. Even if 
one neglects the energy dependence, it is virtually impossible to measure 9 in j directly, 
although experiments have inferred a wide range of characteristic angles in addtion 
to a range of simulation results. Current research is tending to operate under the 
presumption that FI will have to contend with a scenario where the characteristic fast 
electron divergence half-angle is greater than 30°, and possibly even exceeding 50-60°. 

The conversion efficiency from laser energy to fast electron energy, t]l, is not well 
characterized either. A wide range of experimental and theoretical results on this were 
compiled by Davies in (Davies 2009), who noted that the range of results spanned the 
range 10< i]l <90%. Detailed PIC simulations relevant to FI, amongst other results, 
indicate that achieving a tjl in the range of 30-50% is likely, thus making an 'attractive' 
FI scheme still possible provided that efficient coupling to the hot spot can also be 
achieved. 

3.2. Effect of Macroscopic EM Fields 

3.2.1. Return Current and Current Balance As the fast electron beam propagates 
through dense plasma, it will draw a return current that is both spatially coincident 
with the fast electron current density j * and which nearly cancels the fast electron 
current to a good approximation (Bell et al. 2006), i.e. if the return current density is 
jt then, 

j/+jb«0. (6) 

To see how this arises, one can consider the hypothetical case where there is no return 
current. For a wide beam, one can estimate the electric field growth from E ps —jft/eo- 
Since the current densities in FI can easily reach 10 16 A m 2 , one can see that the electric 
field can reach 10 12 V m _1 in 1 fs, which is enough to stop MeV fast electrons on a few 
jum scale. Thus it is clear that a return current will be drawn when the fast electrons 
propagate through dense plasmas. In a fully 3D situation, one might imagine that the 
fast electron current is only globally balanced, but not locally balanced (as in Eq. (6)). 
However, this will lead to the growth of magnetic fields that would destroy the beam, so 
the current neutralization must indeed be co-spatial. The return current phenomenon 
is not particular to fast electron transport in the context of ultra-intense laser-plasma 
physics, and arises in a number of other contexts such as charged particle-beam dynamics 
(Miller 1982) and energetic electron transport in solar flares (van den Oord 1990). 

3.2.2. Resistive Inhibition and Ohmic Heating Current balance implies that j^ ~ — j/, 
and on inserting this into a resistive Ohm's law, one obtains E = — rjjf. The peak 
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resisitivity in many conducting solids will be around 10~ 6 Q m, and even low-Z plasmas 
at a temperature of a few hundred eV will have resistivities of 10~ 7 -10~ 8 Q m. This 
means that the resistively generated electric field can be 10 8 -10 10 V m -1 , which is 
sufficient to inhibit fast electron transport significantly. 

The drawing of the return current also heats the background plasma via Ohmic 
heating with power density j&.E ~ r)jf. From the aforementioned typical values of 
resistivity and fast electron current density this means that the Ohmic heating can 
heat a solid-density target at a rate of 0.1-1 keV/ps. Therefore at solid density the 
Ohmic heating must be included in the energy equation of the background plasma, as 
the heating and thus the effect on resistivity is strong. However at very high density 
(e.g. DT fuel above 100 g cm -3 ) this heating is very small, and thus Ohmic heating will 
not make any significant contribution to the generation of the hot spot. 

3.2.3. Resistive Magnetic Field Generation Current balance also has implications for 
the generation of magnetic field (Bell et al. 1997). An improved estimate for the resistive 
electric field is, 

E = -V3f 4vxB. (7) 

Inserting this into Faraday's law <9jB = —V x E yields 

-XT = ??V x j f + Vr] x j/ + — V 2 B- Vr] x (V x B). (8) 

<jt /J,Q 

The last two terms correspond to resistive diffusion and resistive advection of magnetic 
field, and these are the normal terms that are found in the resistive MHD description 
of a static plasma. The first two terms, on the other hand, correspond to resistive 
generation of magnetic field, and these are due to the presence of the fast electrons. 
Davies noted that one can describe the first term as growing magnetic field which 
pushes fast electrons into regions of higher current density, whereas the second term 
grows magnetic field which pushes fast electrons into regions of higher resistivity. These 
magnetic field growth rates are significant — the magnitude of the growth rate is roughly 
B ~ rjjf/R (where R is the fast electron beam radius). Taking some typical figures 
(jf =10 16 A m~ 2 , i] =10 -7 f2 m, R =5 /jm), yields a growth rate of 2xl0 14 T s" 1 , i.e. 
200 T in 1 ps. Magnetic fields on the order of 100-1000 T will have a significant effect 
on multi-MeV electrons if the fields extend over several microns, insofar as these fields 
can pinch or filament the beam. 

3.2.4- Self- Pinching of the Fast Electron Beam The magnetic field generated by the 
77V x j/ term in Eq. (8) grows in the sense which acts to pinch the fast electron beam 
(Davies 1999, Tatarakis et al. 1998). The tendency of the beam to self-pinch can, 
in principle, be highly beneficial to FI. Counteracting against this self-pinching is the 
angular divergence of the fast electron beam. Bell and Kingham derived a condition 
for self-pinching or collimation for a plasma with Spitzer resistivity by noting that the 
self-pinching condition is R/r g > Q\i 2 with r g the electron gyroradius, i.e. the magnetic 
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field can deflect a fast electron through the characteristic fast electron divergence half- 
angle, #1/2, in the same distance that it takes the beam radius to double. In the limit 
of strong heating (T& S> Tb y i n u), the Bell-Kingham condition (Bell & Kingham 2003) for 
self-pinching is T > 1 where, 

T = 0.13n^ 2/5 (log A)^P^%T°(2 + T 5n )-^RffilL0; a 2 d , 0) 

with n23 being the background electron density in units of 10 23 cm~ 3 , T511 being the 
fast electron energy in units of the electron rest mass, R^ m being the beam radius in 
microns, t psec is the fast electron pulse duration in ps, and 6 ra d is #1/2 in radians. Eq. 
9 shows that the self-pinching is most strongly dependent on the divergence angle of 
the fast electrons, with most other parameters exhibiting much weaker dependence. For 
conditions relevant to FI, self-pinching is marginal and strongly dependent on Q\ii- 

3.2.5. Beam Hollowing Even in a homogeneous plasma, the Vr/ xj/ term can still have 
a significant effect. This occurs in the regime of strong heating as this will produce a 
significant V77 as the Ohmic heating in the centre of the beam is much stronger than at 
the periphery of the beam (oc rjjj). This can lead to the sign of dE x /dr reversing (x| |jf), 
which leads to the generation of a de-collimating magnetic field in the beam centre. In 
turn this will lead to the expulsion of fast electrons from the centre of the beam, and 
this effect is therefore referred to as beam hollowing. Davies first identified this effect in 
(Davies 2003), where he analyzed heating and magnetic field generation in the case of 
a rigid beam model, and he considered different possible resistivity models via i] ex T fe a . 
When a < 1 beam hollowing will eventually occur, as all materials become Spitzer-like 
(a = —3/2) at sufficiently high temperature. 

3. 3. Drag and Scattering of Individual Fast Electrons 

Here we will consider the transport of individual fast electrons through plasmas and 
solids, in other words, we will not consider collective effects arising from the presence of 
more than one fast electron. 

Fast, in this context, refers to an electron traveling at speeds much greater than that 
of the electrons in the material. In this case, the principal effects on the fast electron are 
energy loss and angular scattering. We will present expressions for the rate of energy 
loss, or drag, and the rate of angular scattering and briefly outline their derivations and 
their implications for fast ignition. 

This single particle model will be an adequate description of drag and scattering 
provided that the fast electron density is much less than the electron density of the 
material. To determine exactly how much less requires an accurate calculation of 
collective effects, which we do not have. Work along these lines for the correlated 
stopping of N fast electrons has been presented in (Deutsch & Fromy 1999, Bret & 
Deutsch 2008). In the case of a plasma, this effect should be negligible if the separation 
between fast electrons is greater than the screening length for the fast electron wake. 
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This distance is the dynamical screening length v/u p for v ^> (kT e /m e ) 1 ^ 2 , and the 



plasmas Debye length Ad = JeokT/n e e 2 in the opposite limit. However, there could 
still be significant electromagnetic fields generated by the collective response of the 
material to the fast electrons as a whole that can then be considered independently 
from the drag and scattering. These effects, such as beam-plasma instabilities, are 
discussed elsewhere in this article. 

We briefly note that drag and scattering have received much attention since the 
proposal of FI in 1994 (Atzeni et al. 2009). However, calculations of drag date back 
to the 1930s, with the definitive reformulations of the basic theories being published in 
the 1950s (Sternheimer 1952, Fano 1956). These frequently consider the more general 
problem of drag in matter with bound electrons. Free electrons (namely in conductors) 
are considered in these calculations, therefore they do apply to plasma. These results are 
embodied in (International Commission on Radiation Units and Measurements 1984), 
were recently summarized in (Solodov & Betti 2008, Atzeni et al. 2009), and are 
presented here. The latter two references differ slightly in their angular scattering 
formulas, and in some details of their logic. For a fully quantum-mechanical (but non- 
relativistic) treatment of drag due only to free electrons, including both binary collisions 
and interaction with the plasma medium (e.g. plasmon excitation), see (Ferrell 1956). 

3.3.1. Drag The standard expression for the drag on a fast electron in all matter (solid, 
liquid, gas or plasma, conductor or insulator) is (International Commission on Radiation 
Units and Measurements 1984) 

f = - T^A., (10) 

at A7reQm e v 

T . pv In 2 9 (1/2) ln2 + 1/16 ln2 + l/8 

V7 + Ihujp 2 16 Y 7 

where E, p, and v are the kinetic energy, momentum and velocity of the fast electron, 
respectively. We have introduced the dimensionless parameter L^, which we call the 
drag number. In conventional plasma physics notation it would be called "In A". 
Bremsstrahlung has been neglected. As would be expected, fast electron drag does 
not depend on the velocity or the binding energy of the electrons in the material, since 
we are considering the limit in which these become negligible, it depends only on their 
density n e , which here refers to total, not free, electron density and it is also total density 
that determines the plasma frequency u p in Eq. (11). The value of hu p /e at the typical 
solid density of 6 x 10 28 atoms m~ 3 is Q.lyZ eV, where in this section Z represents 
nuclear, not ionic, charge. In DT, hu p /e = 180^/pioo e V, where p 100 = p/100 g cm~ 3 . 
This result has been extensively tested in cold matter, but not so extensively in plasma 
and never at the densities required for fast ignition, however, there is no reason to believe 
that this lies in a fundamentally different physical regime; drag due to degenerate, free 
electrons is present in metals. 

What changes between materials when applying Eq. (11) is the implication of fast. 
In plasma, it is sufficient for the fast electron to have a velocity a few times higher 
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than the electron thermal velocity (JkT e /m e ) or the Fermi velocity if the electrons are 
degenerate. This will be true in the corona and in the core of ignition targets for all 
cases of interest. In unionized matter, the fast electron must have an energy much 
greater than any binding energy, so in a cone we will have to consider the effect of 
electron binding. Before we do this, we will outline the derivation of Eq. (11) as given 
in (International Commission on Radiation Units and Measurements 1984). 

Fast electron energy loss W above a cut-off W c is calculated using a binary collision 
model and energy loss below W c is calculated using a model for the collective response 
of the electrons in the material. It is assumed that W c is much less than the fast electron 
energy yet much greater than the energy of any individual electron in the material. The 
cut-off W c at which the two models are patched together cancels in the final result, giving 
some confidence that it is accurate, even though neither model is valid for intermediate 
energy losses, for which no analytical model is available. 

For the binary collision model the M0ller cross section is used, an approximate 
solution to the Dirac equation to order av/c (the first Born approximation), where a 
is the fine structure constant, so it includes relativistic effects and quantum spin and 
exchange effects. Experimentally, deviations from this cross section have only been 
detected in close collisions at energies much higher than those of interest here, when 
radiation becomes important. The target electron is assumed to be stationary, which 
requires its velocity to be much less than that of the fast electron, and any binding or 
potential energy is neglected, which requires this to be much less than the energy loss. 
It is also implicitly assumed that the energy loss occurs largely while the electrons are 
close together, because the cross section applies for isolated electrons coming in from 
infinity and being detected at infinity but is being applied to calculate fast electron 
energy loss to only one, immediately adjacent electron among many others. Classically, 
it can be shown that this is an adequate approximation for sufficiently fast electrons 
by considering interaction over a limited distance (Ordonez et al. 1994), and this does 
not represent a significant additional restriction on the theory. We know of no rigorous 
demonstration that this carries over to the quantum case. The calculation follows the 
familiar treatment of binary collisions, with a maximum energy loss of half the fast 
electron energy, since only the fastest electron is followed, giving 



^iw->w, = w^ + ^ - h.2 + —":; t> - _ ^^, (12) 



E 9 (ln2)/2 + l/16 ln2 + l/8 

W c 16 7 2 7 

neglecting terms in W c /E. The first term would be obtained using the Rutherford cross 
section, and the remaining terms represent small quantum corrections due to spin and 
exchange. 

For the model of the collective response of the electrons in the material it is assumed 
that the fast electron moves at constant velocity and that its electric field causes a small 
perturbation of the electrons from their equilibrium positions, so a quantum harmonic 
oscillator model may be used. These are only adequate assumptions far from the fast 
electron and while its velocity changes on a time scale much slower than the collective 
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response time of the electrons. The final result is 



y/2m e v*W c 

L d\ W <W c ~ m fa > l id J 

which when added to Eq. (12) yields Eq. (11). This can be understood in terms of 
energy exchange to plasma waves (plasmons) in quanta of hu p ; that this arises from a 
quantum treatment of electron oscillations in a plasma, which has been considered by 
a number of authors (Ferrell 1956, Bohm & Pines 1953, Pines 1953), is not surprising. 
What is remarkable is that this also arises as the limiting form for fast electrons from a 
general treatment, including electron binding. 

We will now consider the more general case where the binding energies of the 
electrons in the material cannot be neglected, since this will be the case in a cone. 
This leads to the drag being reduced. For a combination of historical and mathematical 
reasons the energy loss due to the collective response of the material is artificially divided 
into two parts and written 



j2p*W c /m e i„2 5 
L d\w<w c = ln j 2^~ 2' ( > 



'- ex 



giving 



-E\ In 2 1 (1/2) ln2 + 9/16 ln2 + l/8 5 



The first part is given by the first two terms of Eq. (14), the basis for which was published 
by Bethe in 1930 (in German). It gives the energy transferred to the excitation of 
electrons by the electric field of a charged particle moving at constant velocity. The 
complexities of dealing with coupled, quantized oscillations of multiple bound electrons 
are hidden in J cx , known as the mean excitation potential, for which there exist a variety 
of theoretical models. In very general terms it can be written as 

lnI ex = J2f lj \n(E J -E i ), (16) 

a weighted sum over all possible transitions of electrons in the material from initial 
energy Ej to final energy Ej, fy being the transition probability. In the simplest 
possible case of a single, undamped, harmonic oscillator of frequency u it is hu. This 
is a good approximation for plasma, giving the mean excitation potential to be hcj p . 
The values normally used for unionized materials are determined by measurements of 
either ion or electron energy loss or of optical absorption, as drag can be treated in 
terms of the absorption of a virtual photon field. Thus it becomes a free parameter used 
to fit experimental data. The reference values are those published in (International 
Commission on Radiation Units and Measurements 1984), available online at (Berger 
et al. 2005). For elements these can be adequately fitted by 9A3Z + 26.1 eV, except for 
hydrogen, where I ex /e is 19.2 eV. For compounds, the stopping due to its constituents 
can be added since chemical structure has been found to have only a small effect on the 
mean excitation potential. The second part is the S, first quantified by Fermi in 1940 
(Fermi 1940) using a purely classical calculation representing the electron response with 
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a single, harmonic oscillator. It gives a reduction in the energy loss due to the electric 
field of the fast electron being shielded by the collective response of the electrons in the 
material, an effect neglected by previous treatments, hence the convention of a negative 
sign (the factor of 2 is another, rather confusing, convention). It is called the density 
effect correction because it increases with electron density. The mathematical reason 
behind this division is the difficulty of giving a straightforward expression for S in the 
general case of multiple bound electrons. It can be obtained analytically in the limit of 
a strongly relativistic electron (Fermi 1940) 



5^2\n 




~7T V — >■ C, 



;i7) 



giving the general result for fast electron drag that we started with. For plasma, where 
J cx = hup, this expression is valid for all cases of interest. For bound electrons, where 
typically I cx ^> hou p , this expression is only greater than for 7 > 1.65I ex /hu p . In 
practice, 7 ^> 1.65I ex /huip is required for Eq. (17) to be a good approximation for 
bound electrons. In solid gold, for example, this requires a fast electron energy much 
greater than 8 MeV. 

Sternheimer has given a simple, approximate formulation of the density effect 
(Sternheimer et al. 1984) that is used in (International Commission on Radiation Units 
and Measurements 1984), which we will write as 

( Eh V 



$ = J2 fn ln 



Ef 



fl B n + fBfn(huj p ) 2 \ \-fh^ 



It 



where f n is the fraction of electrons with binding energy B n (Sternheimer writes this as 
h times the frequency of the absorption edge), /# is 1 if B n = (free electrons) and 2/3 
otherwise, E\ is given implicitly by 

fn{hu p ) 2 1 



£ 



,. flBl + Ef (p/m e c) 2 
and the Sternheimer factor /g is given implicitly by 



J2fn\n\flBl + f B f n (hojp) 



In If 



(19) 



(20) 



which ensures that Eq. (17) is obeyed with experimental values of the mean excitation 
potential, a consistent problem with other formulations. The Sternheimer factor is 
typically between 1.5 and 2.5 (Sternheimer et al. 1984). 

For plasma (B = 0, /# = 1) Eq. (18) and Eq. (19) give the density effect correction 
to be 21n(7) —v 2 /c 2 and Eq. (20) simply gives the mean excitation potential to be hcu p , 
reproducing results we have seen before. 

To illustrate the result for bound electrons let us consider a single binding energy, 
for which it is straightforward to obtain 
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This shows that there is a threshold fast electron energy for the density effect to occur 
in insulators (conduction electrons are treated as free electrons) and that p/m e c should 
exceed fsB/hou p for all electrons before Eq. (11) will be a good approximation, a similar 
constraint to that indicated by Eq. (17). For multiple binding energies a numerical 
solution is required. 

Sternheimer (Sternheimer et al. 1984) gives a 5 parameter fit to the density effect 
correction for numerous elements and compounds. We have found that for Cu and 
Mo the overall drag number from (International Commission on Radiation Units and 
Measurements 1984) is reproduced to within 1% by just using 



*=ln 

2 



1 + -— ^exp(-0.5 

m P c z 1, 



'- ox 



(23) 



which reproduces the limiting forms of the density effect correction, but does not fit 
at intermediate energies. However, here the density effect correction makes a negligible 
contribution to the drag number. This approach could be adapted for insulators with a 
threshold energy E by using 



2 



E — En hu n , n _, 
■exp(— 0.5) 



E > E , (24) 



rn e c i ex 

but we have not verified the accuracy of this approach. 

What is lacking are results for partially ionized matter when the electron is not fast 
enough for Eq. (11) to apply. The only treatment we are aware of is an approximate 
model for the mean excitation potential of bound electrons in an ion, published in a 
difficult to obtain report by More (More 1985), the drag number of the free electrons 
being given by Eq. (11). For the mean excitation potential he used a simplified 
theoretical model known as the local plasma approximation 

In J ex = J f e (r) \n[hu p (f)]dV, (25) 

where f e is the electron probability density function and oo p (r) refers to the plasma 
frequency at the local mean electron density N e f e , N e being the number of electrons. 
In this approximation the mean excitation potential of bound electrons is higher than 
that of free electrons because they are concentrated around the nucleus. To obtain the 
electron distribution around an ion More used the Thomas- Fermi model and found that 
the result could be described by 

T M T ^ exp(1.29(g/Z)°- 72 -°- 18 ^) 

I ex (q) = icx(O) , , (26) 

'l ~ q/Z 



where q is the ionization state. This should be an adequate description for weakly 
ionized many electron atoms when the interaction between electrons of neighbouring 
ions is negligible. It does not give the correct result for hydrogen-like ions, which would 
be expected to have a mean excitation potential of roughly Z 2 times the value for 
hydrogen; More argues that the contribution of this one electron will, in general, be 
negligible. Equation (26) could also be used to estimate the density effect for the bound 
electrons by using Eq. (20) to calculate a new effective value of /s from the new value of 
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J ex , representing an average increase in binding energies, or by using our crude model. 
Binding energies of ions can be measured or calculated, but a limited number of results 
are available and no one appears to have made the effort to apply these to calculating 
fast electron drag. The mean excitation potential and density effect correction will also 
have to be recalculated for compressed material, such as a cone tip in a compressed 
target. The simplest approach would be to start from Eq. (20) to reevaluate the mean 
excitation potential. 

3.3.2. Scattering In solids where scattering is from atoms with a radius a much less 
than the interatomic separation and the de Broglie wavelength of the fast electron is 
much less than a, it is clear that angular scattering can be adequately described in 
terms of binary collisions. An approximate model for the average potential around 
an atom is the familiar exponentially screened potential with a screening distance a 
(Motz et al. 1964, Joachain 1987). Measured and calculated values of atomic radii 
are readily available for all elements; the Thomas-Fermi model gives a simple, general 
result of ao/Z 1 / 3 , where a® is the Bohr radius (5.3 x 10~ n m), although this is not 
accurate for all elements. Using the scattering cross section for an exponentially screened 
potential obtained from the Dirac equation in the first Born approximation (Motz 
et al. 1964, Atzeni et al. 2009), the familiar treatment of binary collisions, integrating 
over all scattering angles, to n, gives 
d(0 2 ) Zn e e A 



dt 2heIp 2 v 



L s (27) 



L s _ a « hi ^ - 0.234 - 0.659^ ^ » 1, (28) 

n c 2 n 

where (8 2 ) is the mean square scattering angle with respect to the electrons 

instantaneous direction of motion, not its original direction of motion, and we have 

introduced the scattering number L s with the -a indicating that it applies to atoms. 

The last term, which is due to the electron spin, had to be evaluated numerically, so 

all terms have been expressed to the same accuracy of 3 significant figures. This differs 

slightly from the expression in (Atzeni et al. 2009) because they calculated (cos 6) not 

(6 2 ). The integral does not diverge at zero scattering angle (infinite impact parameter) 

because in quantum mechanics any potential that falls faster than 1/r has a finite cross 

section for zero scattering angle (Joachain 1987). Here this cross section is approximately 

n[2Zaa/(v /c)] 2 for 2ap/h ^> 1, which surprisingly gives an effective upper impact 

parameter that is smaller than the screening distance a. Since close collisions are not 

modified, using a screened potential in place of the M0ller formula to calculate the drag 

number would make no significant difference. 

The accuracy of the first Born approximation for the exponentially screened 

potential has been carefully analyzed by Joachain (Joachain 1987). He found that 

it is only accurate to order \n(ap/h)Za/(v /c) and only converges for ap/Ti ^> 1 (hence 

our use of this limit), which are quite severe limitations. However, his comparison with 

more accurate solutions shows that what this approximation misses are oscillations in 
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the cross section and that it is accurate for small angle scattering. Since we are only 
interested in the mean scattering angle and the most important factor is the cross section 
for zero scattering, this approximation should not lead to significant errors, with the 
usual provisos that the fast electron energy is high enough for it not to have bound 
states and low enough that radiation is not important. 

In plasmas, following the treatment used for the drag term, we should only use 
binary collisions above some scattering angle 6 C and a statistical treatment of the electric 
field due to random, thermal fluctuations from charge neutrality below 6 C ; hopefully 
C will cancel out. However, there does not exist an adequate model for the effect of 
distant charge fluctuations; all existing models do not deal adequately with interparticle 
correlations due to the electrostatic field and do not include quantum effects, which we 
have seen to be important for distant interactions. The best approach appears to be to 
use Eq. (28) with the Debye length in place of the atomic radius, giving 

L s _ s « In ^-0.234 -0.659^ ^ » 1, (29) 

n c 2 n 

where the -i indicates that it applies to ions. We will now briefly review the theoretical 

models that lead us to this conclusion, in historical order. 

Landau (Boyd & Sanderson 1967) used a series of coupled kinetic equations 

for joint probability densities and set the 3-body joint probability density to zero, 

because it cannot be solved, and obtained an approximate solution for the 2-body 

joint probability density in equilibrium, neglecting particle motion. This showed that 

pairs of particles interact via the exponentially screened potential, but does not prove 

that interactions in a plasma can be reduced to sums over pairs of particles; rather 

it assumes this. Pines and Bohm (Pines & Bohm 1952) used Fourier transforms 

of individual particle positions, but used the random phase approximation, which is 

equivalent to assuming that the particles are uncorrelated. Their treatment went 

beyond that of Landau by considering the effect of particle motion, showing that the 

exponentially screened potential is only accurate for particles with velocities below the 

thermal velocity. Faster particles show reduced, asymmetric screening, for which an 

analytic solution cannot be obtained, but numerical solutions have been published by 

a number of authors (Wang et al. 1981, Decyk 1987, Ellis et al. 2011). Several authors 

have used the Holtsmark distribution for the distant interactions (Chandrasekhar 

et al. 1943, Gasiorowicz et al. 1956), which describes the electric field due to a completely 

random distribution of stationary point charges. This diverges, so an upper cut-off has to 

be introduced. It shows that the net effect of a completely random distribution of charges 

is the same as summing the effect of individual binary collisions with each particle. 

In practice, not all distributions are possible because some will have an electrostatic 

potential energy higher than the total energy of the system. The Debye length, or 

something close to it, then appears as the natural cut-off because it gives the distance 

over which deviations from charge neutrality will give an electrostatic potential energy 

of the order of the thermal energy. Spitzer used two different models (Cohen et al. 1950). 

First, he calculated the random fluctuations in the electric field at a point by applying 
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Poisson statistics to the charged particles in a sphere surrounding it. Like the Holtsmark 
distribution, this assumes a completely random distribution of point charges so diverges 
as the size of the sphere considered tends to infinity and again the Debye length appears 
as the natural order of magnitude for a cut-off to prevent this divergence. He then 
considered the autocorrelation function of the electric field for charged particles moving 
in straight lines, which yet again neglects interparticle correlations and again diverges, 
but this time a cut-off in correlation time is needed. He used 1/oj p , giving much the 
same result as a spatial cut-off at Xd in his first model. This gives a slightly different 
physical picture, with distant interactions being curtailed due to the limited lifetime of 
fluctuations from charge neutrality. A treatment in terms of dipoles has also been tried, 
but not published, and also diverges, although this would not be the case if a quantum 
treatment had been used. 

In summary, these models indicate two practical approaches: 

(i) Sum partial binary collisions over a distance of the order of the Debye length, in 
effect a 1/r potential cut at the Debye length. 

(ii) Sum full binary collisions with all particles using the screened potential. 

We used approach 2 principally because it is more elegant. It also seems reasonable 
to assume that ions move completely at random, which allows us to reduce the many- 
body problem to a sum of binary interactions when considering the mean effect of many 
interactions, because electrons will move to cancel any charge build up before the ions 
are significantly affected by their mutual electrostatic field. The imperfect nature of this 
neutralization due to the thermal motion of the electrons is accounted for by using the 
Debye screened potential, which will apply to the vast majority of the ions since they 
have velocities less than the electron thermal velocity. 

These considerations also lead us to exclude ion shielding in the Debye length, 



which is sometimes included by using Je^kT j{Z + l)n e e 2 in place of JeokT/n e e 2 , a 
conclusion that has been supported by results from numerical modeling (Dimonte & 
Daligault 2008). The Debye length will have to be modified for degenerate electrons. 



A crude approximation is to replace the temperature with JkT 2 + E F where Ep is the 
Fermi energy (Lee & More 1984). 

Unfortunately, it appears that approach 1 will give significantly greater scattering 
than approach 2 because the quantum mechanical result for the screened potential 
effectively cuts off the interaction at a distance significantly less than the Debye length 
(the finite cross section of 7r[2Zaa/(v /c)] 2 for zero scattering). A proper quantum 
treatment of approach 1 is really required, but we can resort to the uncertainty principle 
to iron out this difference; we are considering the interaction of the electron with particles 
within a region of size Ad so they can be attributed a minimum momentum spread of 
order %/2\n, interpreting this as imposing a minimum scattering angle and using the 
small angle approximation gives # m i n ~ h/2p\p,. Using this cut-off with the scattering 
cross section for a 1/r potential obtained from the Dirac equation in the first Born 
approximation (the Mott formula (Mott 1929)) actually leads to a scattering number 
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slightly smaller than approach 2, but the difference is not significant given the crude 
approximation being used. 

We will now consider scattering from electrons, which is normally ignored because 
it is only significant in hydrogen. It is not the same as scattering from ions because the 
maximum energy exchange of half the fast electron energy gives a maximum scattering 



angle of sin -1 a/2/ (7 + 3) and the physical considerations that led us to sum binary 
collisions with all atoms and ions using a screened potential would not appear to apply 
to electrons. For the case of atoms, it seems clear that the electrons do not share the 
same screened potential and that scattering from electrons will only occur while the 
fast electron is inside the atom; the mean effect of the electrons on the total potential 
has been included in the screened potential and we just need to add the effect of the 
irregularities in the potential apparent close to electrons. For the case of a plasma, we 
cannot apply the same argument that electrons are free to move at random and the 
Debye (static) screened potential will not apply to most electrons. The contribution 
of the electrons to the effect of distant charge fluctuations would appear to have been 
included in the screened potential used for the ions, so we just need to include scattering 
due to the random thermal motion of nearby electrons. This amounts to saying that 
approach 1 is more adequate for electrons, with the atomic radius replacing the Debye 
length for atoms. However, we have already argued that both approaches should give 
comparable results, therefore an adequate approximation for electrons should be to 
account only for the reduced maximum scattering angle, giving 

1 1 Q 

L s „ e ~L s - -In 1 —^-. (30) 

The final expression for scattering rate can be written 



d(9 2 ) n e e 



dt 2nelp 2 v 



(Z + l)L s --ln 



1, 7 + 3 



(31) 



2 2 

with L s given by Eq. (28) for unionized material and by Eq. (29) for fully ionized 
material. In (Atzeni et al. 2009) scattering from electrons was dealt with using an 
exponentially screened potential and the result does not differ significantly. 

As with the drag term, we lack results for partially ionized material. In the absence 
of a better treatment, we suggest summing scattering by the ion charge q with a screening 
distance given by the Debye length for the free electrons Xb{q) and scattering by the full 
nuclear charge Z with a screening distance given by the ion radius a(q). This amounts 
to replacing the screening distance in either Eq. (28) or Eq. (29) with \ D (q) q ^ z a(q), so 
as we are only modifying the argument of a logarithm the approximation does not have 
to be particularly good. Ion radii for low values of q are available and for hydrogen-like 
ions it is a /Z, but values for intermediate ionization states are not readily available. 

3. 3. 3. Implications of Drag and Scattering for Fast Ignition We are interested in fast 
electron transport in compressed DT plasma and, for cone-in-shell fast ignition, in the 
cone material, for which gold has been the preferred candidate. When the ignition laser 
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is fired the cone tip will have been heated and shock compressed, so it is not entirely 
accurate to treat it as a cold solid, but we will use these values as an estimate. 

The quantity of principle interest arising from the drag is the stopping distance. 
Assuming the drag number is a constant we can obtain this analytically; 

■ - 4 < E2 (32) 



n e e 4 L d 7 ' 

It is tempting to use the relativistic limit s oc E, but for this to be accurate to within 10% 
requires energies greater than 4.6 MeV, so for most cases of interest the full expression 
should be used. 

Numerical calculations of stopping distances in cold matter are tabulated in 
(International Commission on Radiation Units and Measurements 1984) and are 
available online (Berger et al. 2005). As an example, a 1 MeV electron can penetrate 
up to 400 pm of gold (ps = 0.77 g cm -2 ), so stopping should not be an issue in the 
cone. These calculations also include bremsstrahlung, which allows us to determine 
when this is indeed negligible. In hydrogen energy loss to bremsstrahlung (radiation 
yield) only exceeds 0.1E for E > 100 MeV, while in gold this is reached for E > 2 MeV. 
It is undesirable for an ignition design to have a significant number of electrons above 
1 MeV stopping in a cone, so bremsstrahlung should never be a significant energy loss 
mechanism. 

Fast electron stopping in compressed DT plasma has been considered using Eq. (11) 
in (Atzeni et al. 2009). They give an approximate expression for the stopping distance: 

ps ~ 1.94 ^M£Y O066 cm -2 (33) 

' l + 1.96£ M eV 

which was found to be accurate to within 10% for energies from 1 to 10 MeV and DT 
mass densities from 300 to 1000 g cm -3 . As an example, at 400 g cm -3 a stopping 
distance less than 1.2 g cm -2 requires an energy less than 1.5 MeV. 

Scattering leads to an undesirable increase in the angular spread of electrons, which 
could be quite serious in the tip of a high-Z cone. While energy loss remains negligible, 
the accumulated root mean square scattering angle over a path length s is given by 

02)1/2 « Z#_ [ngj± (34) 

v ' e pv\ 2tt ' v ; 

where n a is atom number density. If we wish to maintain this below, say, 45° (vr/4) for a 

1 MeV electron then for solid gold, tip thickness should be less than 13 /iin (ps =0.025 g 

cm -2 ). Even if the spatial spread of electrons is reduced by a collimating magnetic field 

or vacuum gaps this will not reduce the angular spread, so as soon as the collimating 

effect ends the electrons will diverge. This indicates that a lower Z cone tip would be 

desirable, because even taking into account that thickness should then be increased to 

avoid shock break out roughly as Z" 1 / 2 , the net angular scattering will still vary as Z 3//4 . 

The effect of scattering on ignition requirements for an initially parallel beam of 

electrons entering a uniform sphere of compressed DT plasma has been considered in 

(Atzeni et al. 2009) using a Monte Carlo model. They found that it led to a 10 to 20% 

increase in the energy requirement. 
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3.4- Beam-plasma instabilities 

3.4-1- Motivation Electron beam-plasma instabilities are a long-standing field of 
plasma physics (Davidson 1984). It was early understood that, for a broad parameter 
range, the beam-driven excitation of plasma waves can lead to energy and momentum 
transfer rates between the incident beam and the ambient plasma largely exceeding 
classical (collisional) values (Buneman 1959, Roberts & Berk 1967, Morse & Nielson 
1969, Famberg et al. 1970, O'Neil et al. 1971, Davidson et al. 1972, Lee & Lampe 
1973, Thode 1976, Okada & Niu 1980). Such "anomalous" relaxation or scattering 
processes underlie many scenarios of intense electron beam transport in laboratory 
(Malkin & Fisch 2002, Sentoku et al. 2003) or space (Muschietti 1990, Medvedev & 
Loeb 1999, Achterberg et al. 2007) plasmas. For instance, they were at the basis of 
the pioneering concept of electron beam-driven fusion explored in the 1970s and 1980s 
(Mosher 1975, Lampe & Sprangle 1975, Thode & Sudan 1975, Thode 1976, Miller 1982, 
Sudan 1984, Humphries 1990). Because it relies upon the propagation and dissipation 
of an intense electron current into a large-scale plasma, the fast ignition scheme (FIS) 
has spurred renewed interest in this topic. 

The influence of microscopic beam-plasma instabilities in the FIS could be twofold. 
First, the magnetic turbulence generated by a Weibel-like instability (Weibel 1959, Fried 
1959) in the laser absorption region tends to isotropize the fast electrons through random 
deflections (Adam et al. 2006). As a result, the electrons are injected into the target 
with a large angular spread, which severely constrains the beam energy required for 
ignition: according to Atzeni et al. (Atzeni et al. 2008), the ignition energy increases 
from ~ 25 kJ to ~ 50 kJ when the half-angle divergence of the electron source increases 
from 20° to 40°. Second, the variety of instabilities arising during the beam transport 
could entail an enhanced stopping power which could relax the ignition requirements 
(e.g. (Yabuuchi et al. 2009)). Assuming the beam electrons' mean energy, (£(,), obeys 
the ponderomotive scaling (Wilks et al. 1992), the laser ignition energy, E L , is predicted 
to vary as (Atzeni et al. 2007) 

y300gcm 6 J ^0.5/im r/ L J 

where p is the DT core density, Ao the laser wavelength, r/i the laser-to-electron coupling 
efficiency and fa a parameter (close to unity in the collisional regime) quantifying the 
effective beam range: 

pft = 0.6M£ 6 )gcm- 2 . (36) 

The question therefore arises as to whether the excitation of beam-plasma instabilities 
may entail fa <C 1 so as to significantly decrease El- This could proceed either directly, 
through the unstable wave-beam interaction (Malkin & Fisch 2002), or indirectly, 
through an instability- induced increased plasma resistivity (Sentoku et al. 2003). 

In contrast to past studies, which mostly focused on electrostatic beam-aligned 
instabilities, recent FIS-related theoretical works have considered the whole unstable 
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k-spectrum (Bret et al. 2004, Califano et al. 2006, Bret et al. 2006, Gremillet 
et al. 2007, Bret et al. 2007, Bret et al. 2008, Cottrill et al. 2008, Karmakar 
et al. 2009, Bret, Gremillet & Benisti 2010, Bret, Gremillet & Dieckmann 2010), 
paying particular attention to the quasi- magnetic filamentation modes developing 
normal to the beam direction (Pegoraro et al. 1996, Califano et al. 1998, Sentoku 
et al. 2000, Honda & Meyer-ter-Vehn 2000, Taguchi et al. 2001, Silva et al. 2002, Hill 
et al. 2005, Kato 2005, Tzoufras et al. 2006, Adam et al. 2006, Schaefer-Rolffs 
et al. 2006, Mart'yariov et al. 2008, Polomarov et al. 2008, Karmakar, Kumar, Shvets, 
Polomarov & Pukhov 2008, Shvets et al. 2009, Khudik et al. 2012). Being all the 
stronger when the beam and plasma densities are comparable (Bret et al. 2008), the 
collisionless instabilities are most likely to disrupt the early propagation of the beam 
into the "low"-density regions of the target. Note, however, that from the optimal beam 
intensity found in Ref. (Atzeni 1999), the beam density is expected to be 

^8x 10 21 ( 1 -) (-^-) l cm- 3 . (37) 

VIOOgcm - 3 J \lMeVJ V ; 

Given such extreme values, collisionless instabilities may arise up to solid densities, 
which encompasses the laser absorption region, the cone tip (if any) and part of the DT 
plasma. 

3.4-2. Main instability classes and their related properties Unless otherwise noted, 
we shall restrict our review to uniform, infinite and initially field-free 2-D beam- 
plasma systems. The most general (kinetic) description is afforded by the relativistic 
Vlasov-Maxwell equations, whose linearization yields the following dispersion relation 
for electromagnetic perturbations oc e*( k - x ~ wt ) (Ichimaru 1973) 

(u 2 e zz - k 2 x c 2 ^j (u 2 e xx - k 2 z c 2 ^ - (u 2 e xz + k x k z c 2 ^ = , (38) 

where the dielectric tensor elements read 
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Here, k = (k x , k z ) is the real wave number, u is the complex frequency, u P j = 
(nje 2 /mj€ ) 1 / 2 is the plasma frequency of species j and 7(p) = [1 + (p/rrijc) 2 ] 1 ^ 2 is the 
Lorentz factor. In the following, the index j = (b,p) stands for the electron beam and 
plasma components. Collisional effects are neglected at this stage and will be discussed 
in Sec. 3.4.3. The main ingredient in (39) is the unperturbed distribution function 
fj (p). In the context of the FIS, there is no obvious physical reason supporting a 
particular model distribution for the beam electrons. A variety of descriptions can 
be found in the literature, ranging from monokinetic (Bludman et al. 1960, Pegoraro 
et al. 1996) to Maxwellian-like (Yoon 1989, Yoon 2007, Tautz & Schlickeiser 2005, Tautz 
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Figure 1. Overlay of the normalized growth rate d = ^soj/uj e (shaded colors) and of the 
electric field orientation (arrows) in the (k x ,k z ) space. Maxwell- Jiittner distribution 
functions are considered with nb/n p = 0.1, 75 = 4, T(, = 50 keV, T p = 5 keV and the 
beam drifting along the z axis. 



& Schlickeiser 2006) through waterbag (Yoon & Davidson 1987, Silva et al. 2002, Bret 
et al. 2005, Gremillet et al. 2007, Cottrill et al. 2008) and Kappa (Lazar et al. 2008) 
distributions. However, in order to address potentially large (relativistic) thermal 
spreads, it appears convenient to model the beam-plasma system by means of drifting 
Maxwell- Jiittner distribution functions (Jiittner 1911, Wright & Hadley 1975) 
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rrijC /Tj is the normalized inverse temperature and K2 is a modified Bessel 



function. Two arguments can be made for this model distribution. First, it permits 
an exact resolution of the 2-D fully relativistic spectrum at an affordable numerical cost 
(Bret et al. 2008). Second, it has been shown, under certain conditions, to model with 
some accuracy the relativistic electron phase space observed in laser-plasma simulations 
(Cottrill et al. 2008). Care must be taken, though, in the numerical evaluation of Eqs. 
(39-40) in the complex cu-plane as detailed in Ref. (Bret, Gremillet & Benisti 2010). 

Three instability classes can be identified according to their wave vector's 
orientation and electromagnetic properties. This is exemplified in Fig. 1 which displays 
the k-dependence of the normalized growth rate 

(41) 



8 



3- 



u P 



where u e = [n e e 2 /eom e } 1 ^ 2 is the nonrelativistic total plasma frequency (n e = n^ + n p ) 
for a dilute-beam configuration: nb/n p = 0.1, 7^ = 4, T& = 50 keV and T p = 5 keV. The 
plasma drift velocity follows from the current neutrality condition (3 p = —(3bnb/n p . 

The well-known two-stream modes (Bohm & Gross 1949) are located along the 
beam direction {k x = 0), with a peak growth rate <5 max ~ 0.04 at k Zjraa _ x c/uj e ~ 1/A,. 
These are purely electrostatic plasma waves propagating at the phase velocity Mcu/k ~ 
(3b- Their maximum growth rate is given by the approximate analytical expressions (in 



Fast Electron Transport for FI 



21 



500 " 



200 




0.5 ,1 1.5 

K 0/a> e 






^0 



0.08 
0.06 
0.04 
0.02 



Figure 2. Filamentation growth rate 6 as a function of the tranverse wave vector k x 
and the beam temperature Xt>. Maxwell- Jiittner distribution functions are considered 
with 7b = 2, rib/rip — 0.1 and T p — 5 keV. The cut-off wave vector fcii m [Eq. (44)] is 
plotted in dashed line. 



i-stream 
modes 




Figure 3. (top) Hierarchy of the two-stream, oblique and filamentation modes in the 
(rib/ Tip, 7f,, Tf,) parameter space for Maxwell- Jiittner distribution functions, (bottom) 
Plasma density profiles at the end of the linear phase as predicted by 2-D PIC 
simulations, each ruled by a specific instability class. The plasma temperature is 
T p = 5 keV and the beam flows along the j/-axis. 
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in the hydro dynamic (cold) and kinetic regimes, respectively (Sudan 1984). The 
orientation of the associated electric perturbation can be evaluated from the linear 
relation E x (u,k)/E z (co,k) = (k 2 x c 2 / u 2 e - oj 2 e zz ) / (k x k z c 2 / u 2 e + e xz ) (Bret et al. 2004). As 
expected, Fig. 1 shows that the two-stream modes fulfill k x E = 0. 
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The filamentation instability, which arises in systems composed of counterstreaming 
species, belongs to the family of anisotropy-driven instabilities typified by the Weibel 
instability (Weibel 1959, Fried 1959). Hence, the two designations are often used 
interchangeably in the literature. As the classical Weibel instability, the filamentation 
modes develop preferentially normal to the "hot" (beam) direction k z ~ 0. They 
correspond to aperiodic (JRcu = 0), mostly magnetic fluctuations amplified by the 
repulsive force between electron currents of opposite polarity. In Fig. 1, the filamentation 
growth rate is seen to maximize at k x c/oo e ~ 0.5 with 5 max ~ 0.02. An analytical estimate 
can be derived in the cold limit (T b = T p = 0), which reads (Pegoraro et al. 1996). 



C « hJ— (43) 

V IbUp 

for 1 < k x c/u e < oo. In the kinetic (hot) regime associated to Fig. 1, by contrast, 
the unstable domain is restricted to < k x < k\ im , with the cut-off wave vector (Bret, 
Gremillet & Dieckmann 2010) 



k 2 - 1 - 



r 2 -r a + J(r 2 - T„y + i(F a F 2 - n 
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where T n = J2j(nj/n e )fij/3". Assuming nb/n p <C 1, we have k\ ira ~ 
Pby ( n b/n p )(m e c 2 /Tb). Interestingly, the fastest-growing filamentation wave vector has 
the same scalings as kn m (Bret, Gremillet & Benisti 2010). The shrinking of the 
unstable domain for increasing beam temperatures is illustrated in Fig. 2 for 7, = 2, 
n b/n p = 0.1 and T p = 5 keV. Along with the decrease in /ci im , the peak growth rate 
is found to drop as 5 max oc T b (Bret, Gremillet & Benisti 2010). Further analysis 
shows that, similarly to the cold-fluid scaling (43), the instability is also quenched in 
the high-7fe limit as 5 max oc 7 6 due to the beam's increasing inertia (Bret, Gremillet 
& Benisti 2010). Note that any combination of Maxwell- Jiittner functions with non- 
vanishing /3j's is filamentation unstable (i.e., k\\ ni > 0) due to a finite anisotropy. In 
practice, though, (44) sets an effective stabilization threshold when kn m < 1/L X , where 
L x is the transverse size of the beam-plasma system (typically of the order of the laser 
spot). This incomplete stabilization contrasts with the total suppression occurring 
for model distributions allowing for independent longitudinal and transverse thermal 
spreads (Silva et al. 2002, Bret et al. 2005, Bret, Gremillet & Dieckmann 2010, Tautz 
& Shalchi 2008). In fact, filamentation proves mostly vulnerable to the transverse 
temperature T bx = (p^/7), causing a pressure force counteracting the magnetic pinching 
force. In the simplified waterbag case with weak beam density and temperature, 
stabilization is thus predicted for (Silva et al. 2002) 

fit ~ lbn p ' 
where the beam's transverse velocity spread, A /3b, is related to the transverse 
temperature, T b ± = T bx , through (Silva et al. 2002) 
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2A/3 b n I 1 + A/3, 



(46) 
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which simplifies to T^ x ~ m e '-fbC 2 A(3 b 2 /3 in the limit A/3b <C 1. 

Although the filamentation modes are essentially magnetic, Fig. 1 demonstrates 
that their electric-field component is not purely inductive (k.E ^ 0). This follows from 
the fact that, except for perfectly symmetric systems (i.e, with rib = n p , (3b = —(3 P and 
Tb = T p ), the off-diagonal term e xz in (38) is generally nonzero (Bret et al. 2007). The 
ion response to the resulting space-charge force should therefore be taken into account 
in the weaky-unstable regime (Tzoufras et al. 2006, Ren et al. 2006). 

The spectrum in Fig. 1 turns out to be governed by off-axis modes, thus propagating 
obliquely to the beam. The fastest-growing oblique mode with 5 max = 0.07 is located 
at (k x ,k z ) = (0.8,0.95). As shown in Ref. (Bret, Gremillet & Dieckmann 2010), 
these modes are quasi-electrostatic in a broad system-parameter range including the 
configuration of Fig. 1. For nb/n p <C 1, their maximum growth rate can be estimated 
to be (Rudakov 1971, Sudan 1984, Bret, Gremillet & Dieckmann 2010) 

5° ^KV /3 (47) 

max 2 4 /3 [ lbn J ' W) 

in the hydrodynamic regime defined by 



m P c* 




In the opposite kinetic limit, one has approximately 

Cx-/3^-- (49) 

l b n p 

In both regimes, the longitudinal wave vector of the dominant oblique mode is correlated 
to the dominant two-stream mode (k z c/u e ~ 1/A>), whereas the transverse component 
k x c/u e decreases below unity when moving into the kinetic regime (Bret, Gremillet & 
Benisti 2010). 

The domain of preponderance of each instability class has been computed in 
the (nb/n p ,~fb,Tb) parameter space for a fixed plasma temperature T p = 5keV (Bret 
et al. 2008, Bret, Gremillet & Benisti 2010). The surfaces that delimit regions governed 
by different instability classes are displayed in Fig. 3 and colored according to the 
local maximum (in k-space) growth rate. The two-stream instability prevails for non- 
relativistic beam drift energies (7^ — 1 ^ 1), as well as in weakly relativistic systems 
with hot enough beams. This follows from the quenching of the filamentation and 
oblique instabilities with decreasing $, and increasing T&, respectively. Filamentation 
modes govern systems where the beam and plasma densities are similar (in the FIS, this 
mostly concerns the laser absorption region), whereas oblique modes are dominant for 
dilute relativistic beams. The filamentation-to-oblique transition is mostly determined 
by 75 for dense, cold and weakly relativistic beams, and by rib in the relativistic and 
ultra-relativistic regimes. Note that oblique modes always dominate for hot enough 
relativistic beams. These results are illustrated by the lower panels of Fig. 3 showing 
the plasma density profiles observed in three 2-D PIC simulations, each ruled by a 
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distinct instability class. The spectral characteristics of each modulated pattern have 
been checked to perfectly agree with linear theory. 

3.4-3. Collisional effects Collisions are expected to influence the development of the 
instabilities in the high-density, low-temperature regions penetrated by the electron 
beam, that is, at a distance from the laser absorption region. As a consequence, 
most of the studies performed in this respect have considered dilute collisionless beams 
interacting with dense, nonrelativistic collisional plasmas (Cottrill et al. 2008, Hao 
et al. 2008, Karmakar, Kumar, Shvets, Polomarov & Pukhov 2008, Fiore et al. 2010, Hao 
et al. 2012). Collisional effects are frequently described by simplified Krook-like models, 
which consist in introducing phenomelogical relaxation terms in the Vlasov equation 
(Opher et al. 2002). The most accurate approach of this kind is the particle- number- 
conserving BGK model (Bhatnagar et al. 1954) 



-£ + v.V x / p - e (E + v x B) .V v / P = C(f p 

= -u(5f p -5n p f^/n^) , (50) 

where Sf p — f p — f^ and 5n p = J d 3 v5f p . The BGK model can be generalized to 
conserve momentum and energy as well. A more rigorous approach makes use of the 
Landau collision operators (Shkarofsky et al. 1966, Brantov et al. 2006). In the case of a 
large ion charge (Z ^> 1), electron- ion collisions prevail over electron-electron collisions 
and are described by the operator 

CM,) = ^4^V V i (l - ™) .V./,(») , (51) 



7regmg v \ v 

where / denotes the identity operator and In A is the Coulomb logarithm. In principle, 
the BGK collision frequency v should be adjusted so as to reproduce the plasma 
susceptibility obtained from the Landau operator in the collisional limit, which yields 
v = is e i, where v e i is the usual collision frequency (Bendib 1993). 

An exact evaluation of the collisional two-stream instability using a Maxwell- 
Jiittner-distributed beam and the electron-ion Landau operator has been recently carried 
out (Vermersch & Gremillet 2012). Figure 4 plots the /^-dependence of the growth rate 
for the parameters n e = 10 23 cm~ 3 , n b /n p = 0.01, T b = 100 keV, 7& = 3, T p — 1 
keV, Z = 10 and In A = 2. There follows a collision frequency v e ijbj e = 0.01, that is, 
approximately twice the maximum collisionless growth rate (blue curve). In the presence 
of collisions, the peak growth rate drops from S mgx = 5.3 x 10~ 3 to 1.1 x 10~ 3 , while the 
dominant wave number only slightly decreases. Note that the BGK model (red) yields 
a peak growth rate about 30% lower than the Landau value (black). 

If strong enough, collisions are able to completely stabilize the two-stream modes 
(Singhaus 1964, Hao et al. 2012). This is illustrated in Fig. 5, where the maximum 
growth rate is plotted as a function of the plasma density, the beam density being fixed 
at n b = 10 21 cm -3 . The other parameters are % = 2, T b = 100 keV, T p = 300 eV, 
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Figure 4. Growth rate of the two-stream instability with (red, black) and without 
(blue) e - i collisions for n e = 10 23 cm~ 3 , n b /n p = 0.01, T b = 100 keV, 7& = 3, T p = 1 
keV, Z = 10 and In A = 2. The BGK curve (red) is found to understimate the exact 
Landau curve (black). 
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Figure 5. Maximum growth rate of the two-stream instability for n& = 10 21 cm 3 , 
T b = 100 keV, 7 h = 2, T p = 300 eV, Z = 10 and In A = 2. 



Z = 10 and In A = 2. The exact collisional curve (black) is fairly approximated by the 
expression 



5 w <5 NC - i/ /2 

"max ~ "max ^ei/ ^ 



(52) 

where 5^ x is the maximum collisionless growth rate (blue curve). The relative error 
between (52) and the exact values is found to increase as the instability weakens. 
Complete stabilization (5 max < 0) is achieved here for n p > 3.2 x 10 22 cm -3 , which 
corresponds to v ei juj e > 0.04. By contrast, (52) yields a somewhat underestimated 
stabilization threshold (n p > 2.5 x 10 22 cm~ 3 ). 

Because of the close connection between two-stream modes and oblique modes in 
a broad parameter range (Bret, Gremillet & Benisti 2010), the latter are affected by 
collisions in a similar fashion, exhibiting, in particular, complete stabilization in the 
strong collisional limit (Hao et a!. 2012). 

As first demonstrated by Molvig (Molvig 1975) and further investigated in Refs. 
(Hao et al. 2008, Karmakar, Kumar, Shvets, Polomarov & Pukhov 2008, Fiore 
et al. 2010), a starkly different scenario takes place for the filamentation instability. 
The main result is that for dilute and energetic enough beams, collisions keep the 
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Figure 6. Filamentation growth rate as a function of the wave vector in collisionless 
(dashed lines) and collisional (solid lines) configurations for 7b = 5, T p = lOkeV, 
nb/np — 0.1, v/u) e — 0.5 and increasing beam transverse temperatures: T&j_ = 0.5 keV 
(blue), Tb_L = 9keV (magenta) and T&j_ = 34keV (red) (after Fiore et al. (Fiore 
et al. 2010)). 



system unstable regardless of the transverse beam temperature. Moreover, collisions 
shrink the unstable domain towards small fe's. Figure 6, which is extracted from 
Ref. (Fiore et al. 2010), illustrates these effects by comparing the A;- variations of 
the collisionless and collisional filamentation growth rates for waterbag distribution 
functions with 7 b = 5, T p = lOkeV, nb/n p = 0.1. A BGK collision model is 
employed with u/u e = 0.5. As expected, the instability is weakened and confined to 
decreasing wave numbers as the beam transverse temperature is raised. The instability 
is enhanced in the presence of collisions, especially in the large-temperature limit 
(T b ± = 34keV) where, according to Equation (45), it should be stabilized in the 
collisionless regime. PIC simulations confirm the predicted robustness of the collisional 
filamentation and the generation of filamentary structures larger than in the collisionless 
regime (Karmakar, Kumar, Pukhov & Polomarov 2008, Fiore et al. 2010). Note that 
the highly-collisional filamentation instability corresponds to the so-called resistive 
filamentation instability seen in hybrid simulations (Gremillet et al. 2002, Honrubia 
& Meyer-ter-Vehn 2006, Solodov et al. 2008, Solodov et al. 2009), which is derived 
assuming the return current obeys Ohm's law E = r]j p , where t] is the electrical resistivity 
(Humphries 1990, Gremillet et al. 2002). 



3.4-4- Weibel/ filamentation instability in fast electron generation and transport 
Multidimensional PIC simulations of the fast electron generation in overcritical plasmas 
have shown that the filamentation instability plays a major role in the laser-absorption 
region (Pukhov & Meyer-ter-Vehn 1996, Sentoku et al. 2000, Sentoku et al. 2003, Adam 
et al. 2006, Ren et al. 2006, Okada & Ogawa 2007, Debayle, Honrubia, d'Humieres 
& Tikhonchuk 2010 b). This is so because, for a large enough laser spot (^> Ao) and 
normal incidence, the electron acceleration initially takes place within an essentially 
1-D geometry. As a result of this plane-wave approximation, the transverse canonical 
momentum is conserved: p±(x(t),t) — eA±(x(t),t) = p± . As the vector potential 
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Figure 7. Isosurfacc of the magnetic field |B| associated to Weibel-generated 
fluctuations (averaged over a laser cycle) in a 3-D PIC simulation of a 10 20 W cm - 
laser plane wave interacting with a 50n c plasma. The arrow points along the incoming 
laser direction and the magnetic field is normalized to m e w /e. 



Aj_ vanishes over a few plasma skin depths, the fast electrons quickly recover their 
initial (thermal) transverse momentum |pj_o| (<C \Px\) when penetrating into the target. 
There follows an input electron distribution strongly elongated along the longitudinal 
direction, hence prone to the Weibel/filamentation instability. Magnetic fluctuations 
are then spontaneously generated along the target surface, leading to a fragmentation 
of the fast electron profile into small-scale filaments. 

This process is illustrated in Fig. 7 in the case of a 3-D PIC simulation of a 
10 20 W cm -2 laser plane wave impinging onto a 50n c plasma. At t = lOOcu^ 1 (where coq 
is the laser frequency), magnetic fluctuations have grown in the interaction region to an 
amplitude |B| ~ m e u /e with a transverse size A ~ c/u . The underlying physics can be 
understood as follows. Let us assume for simplicity that the plane wave approximation 
initially holds and that the hot electron distribution is cold in the transverse direction, 
fb(p) oc exp[/iwl + pl/(m e c) 2 ]5(p y ). In the high-energy limit (/i < 1), fi is related 
to the mean Lorentz factor as (7) = K 2 (fJ,) / Ki(fi) ~ 2///. Substitution of the above 
distribution into Eqs. 38-39 (noting that the "hot" axis is now the x-axis) readily yields 
the maximum Weibel growth rate 5 max = u p hJ K (fj) / Ki([i) ~ uJ p bV~ /-thi/x, where 
Upb is the hot electron plasma frequency. Because of the vanishing dispersion in the 
transverse momentum, the growth rate saturates to 5 max for wave vectors k y > ^JJL. 
Assuming that the hot electron density is equal to the critical density (cu p b = uo) and 
that the normalized mean electron energy scales as (7) ~ Oq , where a ~ 2/3 — 1 (Wilks 
et al. 1992, Beg et al. 1997, Ping et al. 2008), one gets 5 max ~ Lu J\na /aQ. As a result, 
the growth rate is comparable to the laser frequency for oq > 1. 

The saturated level of the magnetic fluctuations, -B sa t, can be estimated from 
the widely used trapping criterion (Davidson et al. 1972, Yang et al. 1994, Silva 
et al. 2002, Achterberg et al. 2007, Okada & Ogawa 2007) which expresses the fact that 
the exponential growth phase comes to an end when the electron bouncing frequency 
inside a magnetized filament of period 2ir/k y , w B ~ ^o J (^/l)(k y c/uo)(eB SSbt /m e uo), 
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becomes of the order of the growth rate S(k y ). Using the above estimates with 
(I/7) = K (/jl) / Ki(/j,) and k y = /x 1//2 , one finds the saturated magnetic amplitude 
eB s&t /m e ujQ ~ /i" 1 / 2 = Oq . The maximum quiver momentum being p y = rrie^UB/ky, 
the approximate divergence is p y /p x ~ UB/k y c ~ ^\na . Magnetic deflections 
within the self-generated magnetized filaments then rapidly cause the hot electrons to 
acquire a divergence of the order unity. The Weibel/filamentation instability therefore 
appears as the mechanism mainly responsible for the large angular spread seen in 
simulations (Adam et al. 2006, Debayle, Honrubia, d'Humieres & Tikhonchuk 2010&) 
and experiments (Green et al. 2008). 

The assumption of a zero transverse temperature for the hot electrons actually 
holds only a few skin depths away from the laser absorption region. In reality, however, 
the instability develops within the skin layer, where the electron distribution has a 
finite anisotropy, thus yielding a weaker growth rate. This process has been addressed 
in Ref. (Okada & Ogawa 2007) by fitting the simulated hot electron distribution 
to a semi-relativistic, two-temperature Maxwellian (Okada & Niu 1980). Defining 
A = T x /Tj_ — 1 > and applying the same reasoning as above, the saturated magnetic 
field is expected to scale as 

fn~b A 5 / 2 m e co 



B sa± ft* 0.16a 



2a 

M n c (A + 1)3 



< 0.04a- P^. (53) 

3-D PIC simulations performed with a laser amplitude a^ = 3 and plasma densities 
n e — (1 — 2)n c predict a maximum anisotropy A ~ 2 — 10 and a saturated magnetic 
amplitude in reasonable agreement with the above estimate. 

The 2-D PIC results displayed in Figs. 8(a-d) further depict this self-generated 
magnetic scattering effect in the case of a 10 19 W cm -2 (ao = 3) laser pulse injected 
into a 100n c cone-guided target (Baton et al. 2008). The pulse has a 500^ x duration 
and a I6A0 width. A 1 /im scale-length exponential preplasma is added on the inner 
target walls [panel (a)]. A set of typical electron trajectories inside the absorption 
region are plotted in panel (b). Beside being reflected by the laser field in the low- 
density region (x ~ 200c/w ), the fast electrons undergo strong deflections across the 
skin layer (x ~ 220c/cjo) due to magnetic modulations of amplitude B z ~ 1.5m e Uo/e. 
The resulting electron momentum distribution is shown at t — 300u;^ x before the laser 
maximum [panel (c)] and at the time of the laser maximum [panel (d)]. The root-mean- 
square angle of the fast (> 1 MeV) electrons is found to increase during this time interval 
from (tf 2 ) 1 / 2 ~ 34° to (tf 2 ) 1 / 2 ~ 48°. 

In the plane-wave case, the rapid magnetic build-up breaks the invariance along 
the transverse directions, and hence the transverse canonical momentum is no longer 
conserved. The electron acceleration is modified due to the coupling between the laser 
field and the quasistatic magnetic fluctuations. This multidimentional effect causes 
the transverse velocity of the electrons injected into the target to oscillate at the 
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Figure 8. 2-D PIC simulation of the interaction of a 10 19 W cm~ 2 laser pulse with a 
cone-attached target: initial density profile (a); typical electron trajectories within the 
box shown in (a) and map of the quasistatic magnetic field B z (normalized to m e /woe) 
at the time of the on-target laser peak (b); p x — p y electron phase space around the 
absorption region at t = 300^ 1 before (c) and at the time (d) of the laser peak. See 
text and Ref. (Baton et al. 2008) for further details. 



laser frequency. More precisely, it has been found in Ref. (Adam et al. 2006) using 
a quasilinear analysis that the averaged transverse electron velocity behaves as 



px rx rx 

(vy(x,t)) ^7 _2 y^ / dx I dx" / dx'"kck(x'. 
^ Jo Jo Jo 

xsinlk^- 1 £ t dZA y (Z,t)J , 



X 



(54) 



where A y (t,x) is the laser vector potential, Cfc(xi,x 2 ) = (B k (xi)B_ k (x2)) is the spectral 
density of the perturbative Weibel-generated field B(x,y) = J2kBk(x)e lky . The above 
equation shows that (v y ) changes sign with the laser field in accordance with PIC 
simulations (Adam et al. 2006). 

The late-time dynamics of the magnetized filaments has been frequently 
investigated by means of simulations resolving only the plane orthogonal to the beam's 
axis (Lee & Lampe 1973, Honda & Meyer-ter-Vehn 2000, Sakai et al. 2002, Medvedev 
et al. 2005, Dieckmann 2009). In the case where n^ <C n p , the beam electrons are 
strongly pinched by the magnetic field, which expels the plasma electrons from the 
filament's interior. Further magnetic pinching of the beam electrons generates a strong 
space-charge field accelerating the ions in the radial direction (Honda & Meyer-ter- 
Vehn 2000, Sakai et al. 2002). The nonlinear stage is dominated by the merging of 
magnetically-interacting neighboring filaments (due to incomplete current shielding by 
the plasma electrons), leading to increasingly large filaments. While this merging process 
is accompanied by a steadily-declining total beam current, Polomarov et al. (Polomarov 
et al. 2008) have demonstrated that during its earliest phase, the fusion of sub-Alfvenic 
filaments (i.e., carrying current I < I a — ^bfib^rriec/ [i^e ~ 17.057fc/3& kA) entails an 
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Figure 9. Collisional 2-D PIC simulation of the interaction of a 10 20 W cm -2 laser 
pulse with a cone-guided compressed target: (a) quasistatic magnetic field (in MG) 
at 860 fs; (b) x-profiles of the anomalous resistivity, collisional resistivity, density and 
bulk electron temperature at y = 21 /im (after Chrisman et al. (Chrisman et al. 2008)). 



increase in the magnetic energy and, consequently, a decreasing kinetic energy, whereas 
the coalescence of super- Alfvenic filaments (J > I a) occurring at later times gives rise 
to a slowly-decaying magnetic energy. In the case of comparable beam and plasma 
densities, simulations indicate that the typical filament size increases roughly linearly 
with time as a result of successive coalescences (Medvedev et al. 2005, Dieckmann 2009). 
Provisos, however, must be made concerning the practical relevance for the FIS 
of the aforementioned simulations. First, all of them consider beam electrons with 
weak thermal spreads in contradistinction with the momentum distributions of laser- 
accelerated electrons seen in PIC simulations [Figs. 8(c,d)]. Second, by describing 
the electron dynamics in the plane perpendicular to the beam's flow only, they do not 
capture the parallel or oblique unstable modes described in Sec. 3.4.2. As shown in Ref. 
(Silva 2006) through comparisons with 3-D simulations, the overall influence of these 
multidimensional processes on the beam transport is best reproduced by 2-D calculations 
performed in the plane of the beam's flow. Finally, the simulations carried out in 
Refs. (Lee & Lampe 1973, Honda & Meyer-ter-Vehn 2000, Sakai et al. 2002, Medvedev 
et al. 2005, Polomarov et al. 2008, Dieckmann 2009, Shvets et al. 2009, Khudik 
et al. 2012) employ initially uniform beam profiles with periodic boundary conditions, 
thus neglecting the stabilization provided by the dilution of the diverging beam as it 
propagates away from the injection region. It is then no surprise that a somewhat 
different picture emerges from more realistic simulations of the fast electron generation 
and transport. In particular, for laser intensities 10 20-21 W cm -2 , filamentation is found 
to be confined to the vicinity of the laser-irradiated zone (Adam et al. 2006, Ren 
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Figure 10. 1-D PIC simulation of the interaction of a 10 19 W cm~ 2 laser wave with 
a 100n c , 1 keV plasma: (a) x — p x electron phase space of the interaction region at 
t = 176wq~ ; (b) momentum distribution at various times at x = 98c/ uiq and averaged 
over half a laser wavelength. The solid line plots the high-momentum fit Equation 

(55). 



et al. 2006, Tonge et al. 2009, Debayle, Honrubia, d'Humieres & Tikhonchuk 20105). 
While this region remains (weakly) Weibel-unstable in the nonlinear stage due to the 
destabilizing effect of the ion motion, the interior region becomes stable owing to the 
important dilution of the fast electrons (Ren et al. 2006). Importantly, the filament at ion- 
driven rippling of the target surface triggers additional laser heating mechanisms such 
as the Brunei effect (Ren et al. 2006, Baton et al. 2008). Furthermore, the surface ions 
may be accelerated by the laser radiative pressure to velocities high enough to trigger 
the ion-Weibel instability. The magnetic turbulence thus generated may give rise to a 
collisionless shock of astrophysical interest (Fiuza et al. 2012). 

In Ref. (Sentoku et al. 2003), the magnetized beam filaments have been shown to 
act as random scattering sources for the return current electrons, yielding an effective 
electrical resistivity of the order of z/a = u c /Up€ , where u c is the electron cyclotron 
frequency in the average magnetic field amplitude (\B\). Yet, large-scale laser-plasma 
simulations indicate that this "anomalous" effect only arises within a few microns of 
the irradiated surface, where the backgound temperature is high enough to quench 
collisional processes (Chrisman et al. 2008). 



3.4-5. Electrostatic instabilities in fast electron transport Few studies have addressed 
the influence of the electrostatic (two-stream or oblique) instabilities in the FIS context. 
One notable exception is the simulation work of Kemp et al. (Kemp et al. 2006) 
who showed that, in a 1-D geometry and for a laser intensity of 10 19 W cm -2 , two- 
stream kinetic modes govern the energy transfer from hot to thermal electrons in plasma 
densities < 10 23 cm -3 , whereas they prove strongly inhibited by Coulomb collisions at 
higher densities. 

It is well-known that the nonlinear evolution of the two-stream instability depends 
on the monochromatic or broadband character of the unstable spectrum (Sudan 1984). 
The latter case corresponds to the weakly-unstable, kinetic limit and, to first order, 
is amenable to quasilinear theory (Davidson 1972). Through resonant wave-particle 
interaction (i.e., involving waves satisfying to = k.v), the beam distribution tends to 
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Figure 11. Electron x—p x phase space at t = 800w^~ as predicted by a self-consistent 
PIC simulation (a) and the ballistic propagation of the hot electron source (b). Same 
parameters as in Fig. 10. 



flatten down to the plasma thermal velocity. This weak-turbulence problem has been 
tackled in Refs. (Famberg et al. 1970, Rudakov 1971, Sudan 1984) where the plateau 
formation was found to be disturbed by secondary, nonlinear ion-induced scattering and 
parametric processes. If not collisionally suppressed, this kinetic regime seems to prevail 
in the FIS context due to the broadly-spread and monotonically-decreasing momentum 
distribution of the hot electron source. 

These mechanims are illustrated here by 1-D PIC simulations of the interaction 
of a 3 x 10 19 W cm -2 laser pulse with a 100n c plasma (Gremillet et al. 2012). The 
initial temperature is 1 keV and a lAo scale-length exponential preplasma is added on 
the target surface. In order to obtain a quasi-stationary kinetic energy flux into the 
plasma and, therefore, help identify the unstable beam-plasma processes, the ions are 
kept fixed in a first stage. As a result, the instantaneous laser-to-plasma absorption 
rate has an approximately constant value of ~ 12%. Beyond the laser-irradiated surface 
(x ~ 90c/wo), the x — p x electron phase space displayed in Fig. 10(a) exhibits 2cu Q 
high-energy jets (p x /m e c > 4) typical of the J x B acceleration mechanism (Kruer & 
Estabrook 1985). The electron vortices centered on p x /m e c ~ 1 point to the beam-driven 
excitation of a strongly nonlinear wave close to the absorption region (x < 100c/ u>o). 
This wave, however, rapidly damps out due to bulk electron trapping, hence yielding a 
monotonically-decreasing average momentum distribution, as plotted, at various times, 
in Fig. 10(b). The high-momentum part (p x > 0.2m e c) of this distribution carries a 
density rib/n c ~ 0.2 and, to a good approximation, can be fitted to 

f b (p x ) « 1.67 x lOV 567 * + 0.68e- a5773: (55) 

where ^ x = J\ + p^/(m e c) 2 . Because of its decreasing shape, the source distribution 
is locally stable with respect to electrostatic fluctuations, as also observed by Tonge 
et al. (Tonge et al. 2009). Deeper into the target (x > 105c/w ), though, time-of- 
flight differences generate a transient positive gradient destabilizing the hot electron 
distribution. 

The quasilinear relaxation induced by the two-stream modes developing in the 
dfb/dp x > region is evidenced in Fig. 11, where the PIC-simulated electron phase 
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Figure 12. Energy flux density of the forward-going electrons (in units of m e c 3 n c ) as 
a function of x for a laser intensity of 3 x 10 19 W cm -2 (a) and 10 20 W cm~ 2 (b). 



space at t — 8000;^ 1 is compared to that obtained by ballistically evolving the source 
distribution (55). A plateau clearly forms in the gap region separating the hot 
and thermal electrons, with a width increasing with the distance from the injection 
surface. This proves that the advection time r at j v = x(v^l n ~~ c_1 ) ( wnere ^min is the 
minimum velocity of the hot electrons) is much larger than the characteristic growth 
time tts ~ ( n e/ n b){.&-P x /p z x ) UJ p 1 - This is indeed expected in the present case where 
nb/n e ~ 10~ 3 , p x /m e c ~ 2 — 4, Ap x /m e c ~ 1, v min ~ 0.2c, x ~ 1000c/u p , and hence 
T ndv ~ lOOOujp 1 3> r TS ~ 100w p 1 . Note that the plateau formation is sped up at higher 
laser intensities due to increased beam density. The quasilinear equations describing 
the space-time evolution of the averaged beam distribution function and the spectral 
density of the beam-resonant waves can be analytically solved along the lines of Ref. 
(Zaitsev et al. 1974), by assuming instantaneous plateau formation and using Equation 
(55) for the source distribution (Gremillet et al. 2012). In agreement with the simulation 
results, this model predicts that, for a 10 19 W cm -2 laser intensity, a maximum of ~ 6% 
of the beam energy is converted to resonant waves. Owing to the stable distribution 
source, these waves are subsequently reabsorbed by slower electrons arriving at later 
times. Overall, the wave energy is too weak to affect the beam energy flux. This 
is demonstrated in Fig. 12(a), where the spatial profile of the energy flux carried by 
forward-going (p x > 0) electrons is plotted at various times. Energy is seen to propagate 
at a velocity ~ c with negligible dissipation over ~ 800c/o;o- The spatial variations near 
the right-hand edges of the profiles stem from time-of-flight differences. Figure 12(b) 
corresponds to a 10 20 W cm -2 laser intensity: albeit more strongly modulated than in 
panel (a), the energy flux profiles do not reveal significant dissipation either. These 
findings contrast with the fast relaxation found in the 2-D simulation study of Tonge et 
al. (Tonge et al. 2009). The origin of this discrepancy is not as yet clearly understood: 
it may stem from 2-D physical effects or from the artificial collisionality caused by 
insufficient numerical resolution. 

The time evolution of the wave spectrum in the space region 375 < uqx/c < 400 is 
displayed in Fig. 13(a). As slower and slower electrons reach the detection region, waves 
of decreasing phase velocity are emitted, hence the observed spectral broadening towards 
high fc's. The case of mobile aluminum ions with charge Z = 13 and temperature 
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Figure 13. Electric field spectrum \E x (k,t)\ in the space region 375 < ujqx/c < 400. 
Panels (a) and (b) correspond to collisionless cases with immobile (a) and mobile (b) 
ions, while e — i and e — e Coulomb collisions are described in panel (c). The ion 
charge and temperature are Z = 13 and Ti = 0.2 keV, respectively (other parameters 
identical to those of Fig. 10). 



Ti = 0.2 keV is treated in Fig. 13(b). Weaker and shorter-lived electric fluctuations 
are then generated, as a result of a modulational instability which efficiently scatters 
the beam-excited waves outside the beam-resonant region (Thode & Sudan 1975, Fried 
et al. 1976, Mima & Nishikawa 1984). This parametric process can be modeled assuming 
the primary waves behave as an monochromatic pump wave (ux,ki) decaying into an 
ion wave (oo,k) and Langmuir waves (u\ ± u, k\ ± k). The corresponding dispersion 
relation is (Mima & Nishikawa 1984) 



1 + 



u 2 (kX D ) 2 W E (l + x,)Xe 








(56) 



1 
4 n e T e l + Xe + Xi \D- 
u ± u\) 2 — u 2 — 3v 2 (k ± ki) 2 , Xd is the Debye length, Xj is the jth 



ui p for the wave number 
u p . We have checked that these 



where D± = {u ± u\) 2 — u 2 
component susceptibility and We is the wave energy density. In the present case, one 
has Ui/up = 0.98, ^Xp = 0.53 and WE/n e T e ~ 0.13. Numerical resolution of (56) 
then yields a peak modulational growth rate <5 max = 2.5 x lO^' 
^max = 0.13A5 1 and the real frequency u max = 2 x 10~ 4 
values closely reproduce the simulation results. The high-fc secondary waves generated 
by this instability are strongly Landau-damped by the bulk plasma electrons, which, as 
observed in (Kemp et al. 2006), gives rise to suprathermal tails but negligible ion heating 
(not shown). When Coulomb collisions are switched on, Fig. 13(c) shows that the beam- 
plasma instability is strongly weakened. This is expected since, for the parameters 
under consideration (n e = 100n c , T e = IkeV and T« = 0.2 keV), the collision frequency 
(z/ ei ~ 0.03u p ) is comparable to the collisionless two-stream growth rate. The primary 
waves are then too weak to trigger the modulational instability and the beam-to-plasma 
energy transfer essentially proceeds through the resistive electric field. 

In summary, 1-D kinetic simulations indicate that electrostatic instabilities play 
only a minor role in the energy relaxation of fast electrons generated by 10 19 ~ 20 W cm~ 2 
laser pulses into 100n c plasmas due to the decreasing shape of the electron source 
distribution. Generalization of these results to more realistic 2-D geometries, as 
attempted in Refs. (Tonge et al. 2009, Schmitz et al. 2012), requires further investigation. 
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In particular, the influence of the oblique modes remains to be clearly demonstrated in 
a FIS-relevant laser-plasma setup. 

3. 5. Background plasma physics 

This section discusses the physics of the background medium (i.e. excluding fast 
electrons), that is relevant to fast electron transport and fast ignition. This broadly falls 
in the realm of radiation-hydrodynamics, which we will not review in detail. Instead, we 
focus on aspects of special interest to fast electron transport, which are frequently not 
emphasized in traditional rad-hydro models. These in particular are fluid models that 
include fast electrons, incorporate electromagnetic fields, and account for Fermi-Dirac 
(F-D) statistics (namely background electron degeneracy) in transport coefficients like 
electrical resistivity. 

It is extremely productive to separate transport problems into a background 
medium and fast electrons. This requires distinguishing between fast and background 
electrons, which is generally done based on an intermediate electron energy well below 
that of most fast electrons yet well above that of most background electrons. This 
becomes invalid if the background temperature is comparable to typical fast electron 
energies, either because the background is strongly heated or the fast electrons have 
slowed down significantly. We shall assume that the distinction can be validly made. 

Fast electrons and the background interact via collisions and macroscopic 
electromagnetic fields. Fast electron collisions are discussed in detail in Sec. 3.3 of the 
present article. The e/m fields evolve according to the Maxwell equations, which contain 
the charge and current densities, p and J, carried by the fast electrons and background. 
Since the fast electrons are not atomically bound, it is trivial to find their p and J. 
The background can be much more complicated, depending on whether it is neutral 
matter, a conductor, or a partially or fully ionized plasma. We assume the background 
can be described by a fluid model, meaning it is not a collisionless plasma requiring a 
fully kinetic description. A fluid model applies to neutral matter (with appropriate and 
perhaps difficult models for material properties like equation of state), and for plasmas 
that are sufficiently collisional that the background distribution functions are close to 
equilibrium (e.g. Maxwellian or Fermi-Dirac). 

We focus on the background electron momentum and energy equations, which we 
write in a form following "notation II" of Ref. (Epperlein & Haines 1986): 

<9v f 



^ » e 

dt 


(E — Ec — EjvcO ) 




(57) 


E c = 77 • J e - e" 1 /? • VT e , 




(58) 


Etvc = — - v e x B, 

en e 
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The Maxwell equations including fast electrons are 

- = -VxE, (61) 

eo/io^r = V x B-/i (J e + Ji + J/), (62) 

V-E = eeo 1 f-n e -n/ + 2Z i n i J. (63) 

This section uses SI units and expresses temperature in energy units. Subscripts e, / 
refer to (background, fast) electron quantities. The electron momentum and energy 
equations are written in the rest frame of the ions, so that for instance v e in r]J e should 
be replaced by v e — Vi in a frame where the ions move, cy is the specific electron heat 
capacity, which differs from the Maxwellian ideal-gas result due to e.g. Fermi-Dirac 
statistics. Fluid equations of this type go back at least to Braginskii (Braginskii 1965), 
and require departures from collisional equilibrium to be small. This breaks down, for 
example, when E is large enough that a significant portion of the background electrons 
become runaways, or when v e exceeds the ion thermal speed and triggers the ion acoustic 
drift instability. 

We have expressed the forces as equivalent electric fields (Ec,Ejvc), which arise 
from (collisional, collisionless) effects, respectively. The specific Ec and Ejvc given 
above are those currently implemented in the Zuma code (Larson et al. 2010, Strozzi 
et al. 2012), and neglect certain effects. Namely, E^c lacks the advective term v e • Vv e 
and off-diagonal components of the pressure tensor, and Eq neglects collisions of 
fast with background electrons, rj (resistivity), /3 (thermal force), and k (thermal 
conductivity) arise from collisional or other dissipative effects, which in a weakly-coupled 
plasma are mainly electron-ion (e-i) collisions. They are tensors due to magnetic fields 
or anisotropic distributions, and reduce to scalars for B = and isotropic distributions. 

Transport problems frequently consider situations where background electron 
inertia can be neglected. In plasmas this typically applies for time scales much longer 
than the period of Langmuir waves. Dropping dv e /dt from the momentum equation 
gives an algebraic equation for E in terms of other quantities: E = Ec + Ejvc- We 
call this an Ohm's law. Since this approximation gives E, we cannot treat Ampere's 
law (62) as a time evolution equation for E. Generally the displacement current dE/dt 
is dropped from (62), although it may be fruitful to include it. Ampere's law instead 
gives J e and thus v e , which is no longer specified by the inertialess momentum equation. 
We call dropping both background electron inertia and displacement current the Ohmic 
approximation. Langmuir and light waves are excluded by construction. In addition, 
quasi-neutrality is commonly assumed, which entails dropping V • E from Gauss's law 
and is valid on length scales much longer than the Debye length. One should remember 
that these approximations are independent, even though they are often generically called 
"hybrid models." 

The Ohmic and quasi-neutral approximations currently made in Zuma are 

E = E c + E NC , (64) 
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J e = - j / + /i- 1 V x B, (65) 

n e = ^2 z i n i- ( 66 ) 

Zuma itself does not handle ion motion, but has been coupled to the rad-hydro code 
Hydra (Marinak et al. 2001) as detailed in (Strozzi et al. 2012). 

Regardless of whether an Ohmic approximation is made, a key ingredient is 
specifying the transport coefficients and ionization state. These are aspects where 
ideal (fully-ionized, weakly-coupled, non-degenerate) plasma physics is insufficient in 
transport problems. At high density, it is important that the free electrons obey F-D 
statistics, so that their equilibrium distribution is not the classical Maxwellian. The 
exclusion principle causes the electrons to have random momentum with respect to the 
ions even at zero temperature. This becomes significant when T e < the Fermi energy 
Ep = (h 2 /2m e )(3ii 2 n e ) 2 / 3 . For a sense of scale, fully-ionized solid beryllium (p = 1.84 
g cm -3 ) has Ep = 22.7 eV. Below we sometimes combine T e with Ep in a qualitatively 
correct way, although exact expressions involve F-D integrals. 

A general-purpose framework for transport coefficients in dense plasmas is Lee and 
More's model (Lee & More 1984), which connects the plasma and non-plasma (solid, 
liquid, neutral gas) states. Desjarlais (Desjarlais 2001) provides improvements to their 
model, as well as to the Thomas-Fermi ionization model based on the Saha equation 
and particularly relevant near the metal-insulator transition. Zuma currently employs an 
extended Desjarlais model for ionization and transport coefficients. Whether a material 
at room temperature is a conductor or insulator can be important for experiments, and 
difficult for models geared toward plasmas to capture correctly. 

The main result of the Lee-More model is the electron relaxation time r. From 
this follows the various transport coefficients, including off-diagonal components due 
to magnetic fields. We present the Lee-More model as embodied in Zuma. It has 
been extended to include electron-electron (e-e) collisions along the lines of Refs. 
(Braginskii 1965, Epperlein & Haines 1986). We cast our results in terms of the 
resistivity rj, which has direct physical meaning via the electron momentum slowing- 
down rate v m . For no magnetic field (or the component of rj along B), rj = {m e /n e e 2 )v m . 
t and v m are related by \jv m = rA a . A a accounts for electron F-D statistics and involves 
F-D integrals F n : 

A° = 4 F2 ^ m\ 

3[H-exp(-/})]F 1/2 (/})2' l ) 



ot 



(x) = dtt r 
Jo 



e f - x + 1 



-i 



(6* 



A a is given in Eq. (25a) of Ref. (Lee & More 1984), except with a typographical error 
that F 2 incorrectly reads F 3 there, ft = /i/T e where // is the electron chemical potential, 
defined implicitly by F 1/2 (ft) = (2/3)#~ 3/2 with 9 = T e /E F . Antia (Antia 1993) provides 
rational function approximations to F n for several half-integer orders and their inverses, 
which we use to directly find fi(6). 
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A formula for r that spans the plasma, neutral-gas, and condensed regimes is 

r = max(r ec , r me i t , r min ) . (69) 

r" 1 = r" 1 + r en 1 defines the electron collision time off both charged ions (r ei ) and neutral 
atoms (r e „), with rates added. r me i t and r min stem from a Bloch-Griineisen melting model 
(Ziman 1961), and a minimum time based on inter-atom spacing Ri = (3/47rnj) 1 / 3 . 
Tg" 1 = n n a n v where n n is the number density of neutral atoms, a n the cross-section, 
and v = 3 • 2~ 1 / 2 (T e /m e ) 1//2 # 3 / 2 .F 1 (/i) the average electron speed. The limiting values of 
■u 2 are (8/n)T e /m e for 3> 1 and (9/8)eF/m e for <C 1. Approximately, r min ps -Rj/w, 
and Desjarlais has discussed refinements to this (Desjarlais 2001). The melt model gives 
Tmeit ~ 50(T me it/T e )r min , with the material-dependent constant 50 decreasing somewhat 
for T e > T me it (see Lee and More for details). The melt model applies to conductors with 
strong ion correlations, such as a periodic lattice, and not to insulators or gases. In the 
periodic case, the electron wavefunction becomes a Bloch wave in a periodic potential, 
and essentially does not undergo Coulomb collisions off single ions. Instead, electrons 
slow down due to interactions with phonons. As temperature increases and the ions 
become uncorrelated, Coulomb collisions with ions dominate, and r ec applies. Taking 
the maximum of the three r's is a crude way of capturing the real, more complicated 
physics. 

We now discuss r e j, which falls closest to the realm of traditional plasma physics. 
We consider one electron species colliding with one ion species; for multiple ion species 
the collision rates and therefore the 77's add. Lee and More find 



77 1 r e 2 m l J 2 n. 

r] 39 3 F 2 {[i) eJ n e 



//„ = -KV2—^5 ee -Z 2 In A ei . (70) 



l + iAn-^Oy] V2P . (71) 



770 is the fully-degenerate result, and 5 ee accounts for e-e collisions. In the non- 
degenerate (Spitzer) limit i] = i] S with i] S /t]q = (n 1 / 2 /8)9~ 3 / 2 , or r] S oc T^ 3 ^ 2 (neglecting 
temperature dependence of 770 via Z$ or In A ei ). An approximate form with the correct 
small and large 6 limits is 
rj_ 
Vo 

For p = 1.72 the relative error in i]/i] is at most 2% for all 9. Some workers 
approximately include Fermi degeneracy by capping rj at the value at some temperature. 
The Spitzer r/s equals the fully degenerate r/ for 9 = (it 1 / 2 /8) 2 ^ 3 = 0.366. 
We use the Lee-More Coulomb logarithm: 

1 



In A P i = max ( 2, - In 



1 + A Z 



(72) 



with A = b max /b m i n the coupling parameter; A ^> 1 for a weakly-coupled plasma. 
bmax = [^DH + ^f] 1 ^ 2 1S ^ ne overall screening length. The Debye-Hiickel screening length 
is given by 



K Z h = - ^H-E 1 ^ (73) 
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Figure 14. Resistivity for beryllium at density 1.84 g cm -3 from various models. See 
text for details. 

with T = (T 2 + (4/9)-E^) 1,/2 . The minimuni impact parameter is 

h 2 = h 2 4- h 2 (7A\ 

u min min,Q > "min,Ci V* / 

(75) 
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b m in,Q is the de Broglie wavelength, and b m i n ,c is the classical distance of closest 
approach. 

The Lee-More model has proven successful at capturing the results of experiments 
or more detailed models. Most concrete implementations of Lee-More involve several 
material-dependent adjustable parameters. They can frequently be chosen to replicate 
more correct results. Like any semi-analytic model, Lee-More can be applied over 
wide parameter ranges and usually gives smooth results. Tabulated output from more 
detailed models can include more physics. The typical drawbacks of tables include the 
limited parameter range over which they were generated, and the difficulty of tabulating 
a high-dimension domain (for instance, a dopant ion species of variable concentration 
increases the table dimensionality). 

More sophisticated models than Lee-More exist, and are particularly necessary 
in the non-plasma regime. One is the Purgatorio code developed at LLNL (Wilson 
et al. 2006, Sterne et al. 2007). It solves the Dirac equation for bound and continuum 
electron states surrounding a single ion. Transport coefficients like resistivity are found 
using an extended Ziman formulation (Ziman 1961, Evans et al. 1973, Rinker 1988), and 
require specification of the ion correlation function, or equivalently the ion structure 
factor. See Ref. (Hansen et al. 2005) for details. A similar approach to transport 
coefficients is presented in Ref. (Rozsnyai 2008). 

Figure 14 plots rj vs. T e for 1.84 g cm" 3 beryllium from various models. The solid 
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Figure 15. For beryllium: (a) Where different terms dominate b m i„ and b max . (b) 
Coupling parameter A with black curves for A = 2,5, 10,30, 100,500. White dashed 
curve is fit to A = 10 contour: T e = 1.5 keV(p/500g cnr 3 ) - 45 . 



black curve is from Purgatorio f. The ion structure factor used was a combination of 
the results of Baiko et al. (Baiko et al. 1998), the one component plasma model, and 
Debye-Hiickel theory. Purgatorio calculates the charge state as well, which increases 
smoothly with T e , from 1.5 at room temperature, to 3.2 at 100 eV, and asymptotically 
approaching 4 for higher T e . We use Purgatorio's Zj in the other calculations. Since 
Zi > 1 for all T e we neglect electron- neutral collisions (r en — > <x>). The solid red curve 
comes from Lee-More's r e j for e-i collisions, modified to include e-e collisions (5 ee ^ 1). 
The dashed red curve is the Spitzer i] S (E F /T e — ¥ 0). The solid blue curve is the full 
Lee-More model, with numerical parameters chosen to give a decent agreement with 
Purgatorio at low T e , and with e-e collisions neglected (5 ee = 1). This last choice gives a 
slight difference between the red and blue curves at high T e . It is easy to include S ee , but 
we omit it to demonstrate its magnitude. The dashed blue curve comes from just the 
melting model r = T me i t . Although we have found parameters that bring the Lee-More 
model into decent agreement with the more complete Purgatorio results for the chosen 
density, those values are likely not optimal for all densities. 

Some recent attention has been paid to transport coefficients for strongly coupled 
plasmas, i.e. A ~ 1 (Baalrud 2012, Glosli et al. 2008, Daligault & Dimonte 2009). 
These usually do not consider F-D statistics and connection with the non-plasma state, 
so that a strongly-coupled generalization of Lee-More is not yet at hand. Following 
Baalrud (Baalrud 2012), strong coupling significantly affects the Coulomb logarithm 
when A < 10. We assess the potential importance of strong coupling for beryllium, 
using the Purgatorio Zj, in Fig. 3.5. Panel (a) indicates where different terms dominate 
b m i n and bmax-, and panel (b) plots A including a fit to the A = 10 contour. Strong 
coupling is estimated to be significant below this contour, which includes regimes of 
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interest to fast ignition and high energy density physics more generally. 

4. Simulation Methods 

In this section we will review the simulation methods that have been developed for 
studying fast electron transport, and particularly those that are used in the Fast Ignition 
context. 

4-1. Vlasov-Fokker-Planck Codes 

The Vlasov-Fokker-Planck (VFP) equation for electrons describes their motion through 
phase space under the action of the local average Lorentz force and the microscopic 
field fluctuations that give rise to small-angle collisions with other electrons and ions. It 
is usually expressed in Cartesian geometry without giving the details of the collisional 
term (Thomas et al. 2012): 

| +v f- e(E + VX B,.|.(|), 

where x and p are the phase-space position and momentum coordinates (respectively), 
/ = / (x, p, t) is the electron distribution function, E + v x B is the Lorentz force and 
the term on the RHS accounts for the scattering in momentum-space due to collisions. 
When this equation is solved by the use of computational particles it is known as 
"collisional Particle-In-Cell" (collisional-PIC). The collisional-PIC technique has many 
advantages: robustness; good momentum and energy conservation; no stability or 
magnitude restrictions in momentum-space; very accurate advection in momentum- 
space; and it naturally concentrates computational effort in well-populated regions of 
phase-space. The collisional-PIC technique also has the great benefit of over 40 years of 
research experience behind it. However, it does suffer one major drawback in that it also 
introduces noise into the numerical result. The effect of this noise is arguably not yet 
well-explored in the regime of the dense plasmas that arise in fast-ignition research. 
For this reason a small number of codes based on finite-difference techniques (Bell 
et al. 2006, Robinson & Sherlock 2007, Sherlock 2009) have arisen in recent years. While 
techniques based on finite-difference in phase space eliminate noise, they do not currently 
possess most of the aforementioned advantages inherent in the particle approach (and in 
fact are often deleteriously affected by their converses). Numerical diffusion also occurs 
(although modern techniques do minimize it) and this may adversely affect the physics. 
Nevertheless the finite-difference approaches have much merit and also serve as a useful 
"reality-check" on the results gained with the particle techniques. 

Although the VFP equation is largely valid in Fl-relevant plasmas, it should be 
kept in mind that the collision term requires corrections of order 1/lnA (the inverse 
Coulomb logarithm), which should be necessary when the plasma is initially cold and 
dense. Since the VFP equation is valid over all of momentum space, it is possible to 
solve it for the distribution function of the energetic particles only, and this is the basis 
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of the VFP-hybrid technique (as used in e.g. (Robinson & Sherlock 2007)), where the 
background electrons are treated as a simple fluid (see the section on hybrid methods). 
In fact due to the heavy computational demand of solving the VFP equation, the hybrid 
technique is by far the most common approach. 

Since angular-scattering is important for electrons, it is advantageous to use a 
spherical coordinate system in momentum-space, as angular scattering can be easily 
expressed as diffusion of the distribution function in the angular coordinates. Recasting 
(77) in spherical-coordinates in momentum-space (p, 9, 4>) and introducing diffusive and 
drag-like (i.e. in p) collision terms gives: 

— — h v cos v sin (p— — h v sin 9 sin (p-^- 
at 




where 



+\ (Y ee n e + Yim) ^- [ -—r-~ + -—— { sin <>~: } \ - (78) 



F x = —e (E x — vB z sin 9 sin 0) , (79) 

F y = -e (E y + vB z cos9 sin (f)) (80) 



are the components of the Lorentz force and v = p/'jm, 7 = J\ + p 2 /m 2 c 2 , Y ee = 

Aw (e 2 /4-7reo) In A ee and Yi = 4tt (Ze 2 /47reo) In A e j. The collision terms on the RHS are 
taken from (Johzaki, Nakao & Mima 2009) and are valid for hybrid VFP simulations 
only: the first term accounts for the dynamic friction of fast electrons with the cold 
electrons and ions, while the second term accounts for the angular scattering of fast 
electrons off cold electrons and ions. 

Equation (78) can be readily solved in flux-conservative form, as is done in the 
FIDO simulation code (Sherlock 2009), which uses Piecewise- Parabolic-Interpolation to 
compute the fluxes combined with the Finite-Difference Time-Domain (FDTD) scheme 
for Maxwell's Equations to compute the fields (see e.g. (Villasenor & Buneman 1992)). 
This treatment is particularly advantageous when the field (acceleration) terms 
dominate, as is the case when absorption in strong laser fields is modelled. 

An alternative form of Equation (78) is possible by expanding the distribution 
function in momentum-space in a spherical-harmonic basis: 

N n 

/(t,x,v) = £ £ f?(t,x,v)P? (cos 9) e im * (81) 

n=0 m=—n 

where 9 is the angle between the velocity vector and the spatial coordinate and the P™ 
are the associated Legendre functions (Bell et al. 2006). This gives rise to a large set of 
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coupled partial differential equations for the coefficients: 

bq_ + Mn+v i {Y ^ + Yinj) ^ _ Cm = (82) 
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fl^i 1 + (n + m + l)(n + m + 2)flK 1 ]}, (87) 

m lin + o L J J 

where G™ = df™/dp — nf™/p , if™ = df™/dp + (n + 1) /™/p and C ee accounts for 

the friction between fast electrons and the background plasma. Equation 87 is for the 

ID case only (see (Bell et al. 2006) for the full 2D equations). This form of the VFP 

equation has the advantage that it can be solved with relatively fast and simple (for 

example Runge-Kutta) integration schemes, provided the driving fields are small (in 

comparison to the laser fields). Numerical schemes that can handle large perturbations 

to the distribution function are complex and slow. It also allows for a particularly 

accurate treatment of the magnetic field terms, which are reduced to algebraic form when 

differenced. This form for the VFP equation was used in, e.g. (Bell & Kingham 2003). 

4-2. Hybrid Ohmic Codes 

It is apparent from the disparity between the cold and fast electron characteristics, as 
described in Section 3 that the problem of fast electron transport is computationally 
'stiff' (disparate length and time scales), and that this also allows a natural separation 
of the problem into two interlinked models. It is this observation that has motivated 
the development of the 'hybrid' code. The term 'hybrid code' appears in many places 
in plasma physics, and the term often denotes very different things. In the case of 
fast electron transport, the term denotes a code in which a kinetic treatment is applied 
to a distinct fast electron population, and a fluid treatment is applied to a distinct 
background plasma. It is frequently assumed that the background plasma will respond 
instantaneously to the fast electrons to ensure quasineutrality, which is valid on length 
scales larger than the Debye (or other screening) length. An independent assumption, 
which we call the Ohmic approximation, is that the electric field can be determined 
from a suitable generalized Ohm's Law, with displacement current c^E dropped from 
Ampere's law. In all cases the magnetic field is evolved from Faraday's Law. 

This splitting of the populations will only be a good approximation when n/ <C n&, 
and when the fast electron mean energy is very much greater than the mean thermal 
energy of the background electrons. In the case of the ultra-intense laser-generated fast 
electron transport problem these conditions will be quite easily satisfied at material 
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densities above 1 g cm~ 3 . Although the fast electron population is very much less 
dense than the background electrons, the fast electron current density is still sufficiently 
large to generate electric and magnetic fields with energy densities comparable to the 
fast electron energy density unless there is a compensating return current carried by 
the background electrons (see §3.2). This leads to the conclusion that, to a good 
approximation, jj + j b ?» 0. A more accurate approximation is j/ + j& ~ A^V x B. The 
key equations for E and B in the hybrid approximation then become, 

E = -m f + — Vx B, (88) 

and 

— - = V V x j f + Vr? x u + — V 2 B Vr? x B. (89) 

at /i /i 

The Ohm's Law shown in Equation (88) can be easily extended to include a number of 

additional terms (see Sec. 3.5 of this review) — here it is just given in one of the simplest 

forms. The kinetic treatment of the fast electrons is mathematically described by a 

suitable kinetic equation, i. e. (77). Usually this is solved by using the standard Particle- 

in-Cell methods, but with collision operators for the angular scattering from background 

ions and electrons, and drag from background electrons, included via a Monte Carlo 

method. The background plasma is described, in general, by a set of hydrodynamic 

equations, although for some problems it is reasonable to treat the background plasma 

as essentially static. Even if the background plasma is static, its temperature must 

evolve due to both Ohmic heating of the background electrons, and collisional drag on 

the fast electrons. These effects must be incorporated into the background electron 

energy equation. 

This must also be accompanied by a prescription for the resistivity. This can be 
from a theoretical model (Spitzer resistivity or Lee- More), an empirical model, or even a 
heuristic model. Although the background resistivity is not calculated self-consistently 
(unlike a purely kinetic model), one advantage of the hybrid approach is that it is 
relatively easy to use a resistivity model that better treats the 'warm dense matter' 
regime. This regime is unavoidable both in solid density interactions and Fast Ignition, 
as the resistivity will only be very well approximated by the Spitzer resistivity well above 
the Fermi energy. For DT at 1 g cm~ 3 this would be temperatures above 14 eV. 

Although the hybrid approach makes a number of approximations it also has a 
number of powerful advantages. Firstly the fluid treatment of the background means 
that the very small length and time scales of the background plasma can be ignored, 
and thus much larger time-steps can be used which allows large problems to be run 
quickly. Secondly the model is very computationally robust. Thirdly hybrid codes 
are easy to write and maintain. Fourthly, hybrid codes allow a lot of physics to be 
included easily. There are a number of such hybrid codes both in current use and 
reported in the literature. This includes the unnamed code of Davies (Davies 2002), 
PETRA(Honrubia & Meyer-ter-Vehn 2006), LED A (Robinson & Sherlock 2007), and 
ZEPHYROS (Robinson et al. 2012). 
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4-3. Hybrid Implicit PIC Codes 

LSP (Large Scale Plasma) is a one-, two-, and three-dimensional particle-in-cell (PIC) 
code developed by Mission Research Corporation (Welch et al. 2001) and currently 
maintained by Voss Scientific (Welch et al. 2006, Thoma et al. 2011). LSP has a 
hybrid mode, as well as several types of electromagnetic field solution available: the 
standard (explicit) leapfrog algorithm and implicit algorithms using iterative Alternate- 
Direction Implicit (ADI), a two-step ADI, and matrix inversion. The implicit algorithms, 
particularly the two-step ADI, are useful in relaxing the Courant limit on the time step. 
An iterative electrostatic algorithm is also available for simulations in which fields are 
slowly varying. For short-pulse LPI and fast electron transport problems, the relativistic 
implicit-PIC code ELIXIRS (Drouin et al. 2010) employs a similar approach. 

LSP has several options for advancing particles: the standard momentum- 
conserving and energy-conserving PIC algorithms, cloud-in-cell (CIC) algorithm, and 
direct implicit particle/field algorithm (Hewett & Langdon 1987) which can be used in 
either the PIC or CIC models. The direct implicit algorithm is used most often. The 
benefits of the direct implicit algorithm are that the usual charged particle limitations 
on the time step, namely the need to resolve the cyclotron and plasma frequencies, 
are relaxed although both frequencies cannot be under resolved at the same time and 
position. The implicit algorithm is useful for very dense plasmas (that occur in fast 
ignition studies) where the details of electron plasma oscillations can be ignored. By 
damping unresolved high-frequency Langmuir modes the direct implicit algorithm allows 
stable modeling of dense plasmas without needing to resolve these fast modes. 

The energy conservation as well as the speed of the direct implicit calculation is 
further improved by including a nonrelativistic inertial fluid model for the electrons 
in which the directed and thermal energy of the electrons are treated separately. The 
equation of motion for the fluid electrons is of standard Braginskii type (Braginskii 1965). 
It includes a frictional force to model collisions with other particle species. For 
temperature, the new energy equation for an ideal gas is added including the pdV 
work, energy exchange between species, thermal conduction, and Ohmic heating rate. 
Inelastic losses with neutrals can also be included. In some circumstances, kinetic effects 
become important. Examples include runaway, where a hot electron population coexists 
with a thermal one, or acceleration of less dense electrons or ions from a biased plasma. 
In the hybrid mode, LSP permits dynamic reallocation of particles between the fluid 
and kinetic description. Fluid particles with directed energy much greater than thermal 
energy transition to kinetic particles. Kinetic particles with energy less than the ambient 
fluid thermal energy transition to fluid particles. These transitions result in energetic 
electrons treated kinetically and dense thermal plasma electrons as a fluid. 

The LSP simulations include an algorithm to model electron-electron, ion-ion and 
electron-ion collisions. For kinetic particles, this involves first constructing drifting 
Maxwellian distributions at each grid cell. A particle of a given species is then 
elastically scattered isotropically in the center-of-mass frame off a particle obtained 



Fast Electron Transport for FI 46 

by sampling this distribution. Collisions between different species (both kinetic and 
fluid) are separated into an energy push and a frictional momentum push. The energy 
and momentum transfers from one species to another are accomplished by summing the 
changes from each interaction on the grid. The collision frequencies are determined from 
the Spitzer formulation. Optionally, the more accurate Lee-More (Lee & More 1984) 
model with Desjarlais corrections (Desjarlais 2001) is available for collisions of fluid 
electrons and background ions with the ion charge state calculated with a Thomas- Fermi 
equation-of-state with pressure ionization corrections. Monte-Carlo type scattering 
model with the drag and scattering formulas of (Solodov & Betti 2008, Atzeni et al. 2009) 
is also available for collisions of kinetic electrons and background plasma electrons and 
ions. 

The LSP hybrid implicit approach to dense plasma modeling is alternative to 
that of hybrid Ohmic codes, described in the previous subsection. In those codes, the 
background plasma is modeled as a collisional fluid and charge neutrality is assumed. 
The electric field is found from Ohm's law and the background return current is 
found from Ampere's law without displacement current. This reduced-model approach 
is inapplicable to laser-plasma interactions, or low-density regions with, e.g., Debye 
sheaths. The LSP model solves full Maxwell equations with displacement current and 
is applicable in the laser-plasma interaction and low-density regions, provided a kinetic 
description for plasma electrons is used there (while a fluid description can be used 
elsewhere in the same run). Time steps and cell sizes can be chosen that resolve the 
laser-plasma interaction (LPI) near the critical surface. These time steps, while explicit 
in the LPI region (electron density, n e = 1.1 x 10 21 cm' 3 for 1 /im laser light), can be still 
highly implicit to the plasma frequency in the solid density target (n e >10 23 cm" 3 ). The 
entire process of laser propagation into the underdense plasma, fast electron production 
and transport into the dense plasma, and ion acceleration from the plasma-vacuum 
interface can be modeled with LSP (Welch et al. 2006). 

Detailed equation-of-state (EOS) and multi-group diffusion radiation transport 
modeling capability was recently implemented in the LSP framework with fluid particles 
(Thoma et al. 2011). The EOS and opacity data needed for the algorithm are pre- 
calculated by the Propaceos code (Macfarlane et al. 2006), which utilizes detailed atomic 
models for plasmas in local thermodynamic equilibrium (LTE) as well as non-LTE states. 
The EOS model is used to evolve the ion charge state and introduce non-ideal gas 
behavior. The radiation energy density field is calculated, which is coupled to the 
plasma. 

5. Review of Ignition Scale Calculations 

This section summarizes the FI calculations carried out so far, showing the dependence 
of the ignition energies on the electron beam parameters. Electron-driven FI modelling 
relies on the characterization of the fast electron source, the transport from the 
generation zone to the compressed core and the energy deposition in the fuel. Thus, 
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the complete description of FI requires the integration of different models/codes that 
deals with different spatial and temporal scales. Fully integrated calculations are not 
possible with the present computer resources. Here, we focus our attention on the 
partially integrated calculations that consist of characterizing the fast electron source via 
experiments or PIC simulations and using this source to perform fast electron transport 
calculations coupled to radiation-hydrodynamics, including fusion reactions. This 
'integrated' model has been used so far to estimate the electron beam requirements in 
the fast ignition scenario (Solodov et al. 2008, Honrubia & Meyer-ter-Vehn 2009, Strozzi 
et al. 2012). Relativistic Fokker-Planck models for electron transport in sub-ignition 
targets have been developed also within the context of the FIREX-I project (Yokota 
et al. 2006). 

5.1. Ignition energies of perfectly collimated electron beams 

We assume for the moment that a perfectly collimated beam impinges on a DT assembly 
at a time close to the peak compression. Here, we do not take into account the fast 
electron scattering nor the EM fields generated by the electron beam. We first consider 
the target proposed by Solodov et al. with a DT mass of 0.5 mg compressed by a 300 kJ 
nanosecond laser pulse to a peak density of 500 gcm~ 3 (Solodov et al. 2007). Assuming a 
mono-energetic electron beam of 20 /im radius and a pulse duration of 10 ps, the lowest 
ignition energy E ig = 16.2 kJ is reached for 2 MeV electrons (Solodov et al. 2007). This 
ignition energy is much higher than the 7 kJ obtained from Eq.(l) for p= 500 g/cm 3 due 
to the target density profile and also because the beam radius and pulse duration do 
not have the optimal values given by Eqs.(2-4). For the more realistic case of electrons 
with a relativistic Maxwellian energy spectrum, the lowest ignition energy raises to 
21.5 kJ for an electron mean energy (E) = 1.25 MeV. This energy is substantially 
lower than the 2 MeV found for mono-energetic beams because Maxwellian electrons 
deposit their energy over a larger region. In addition, it is much lower than those 
obtained in experiments and PIC simulations for laser intensities around 10 20 Wcm~ 2 . 
More realistic electron energies can be obtained by using the ponderomotive scaling 
formula (E)/m e c 2 = [1 + JlA 2 /1.38 x 10 18 ] — 1 , which relates the laser intensity 
II (in Wcrn" 2 ) and wavelength A (in /im) with the electron mean energy (E) (Wilks 
et al. 1992). For instance, let us consider a Gaussian laser pulse with a duration of 
13.8 ps and A = 1.054 /im impinging on the fuel assembly mentioned above. Assuming 
that the electron mean energy is given by the ponderomotive scaling and a laser-to-fast 
electron conversion efficiency of 50%, one obtains (E) = 6.3 MeV and an ignition energy 
E ig = 53 kJ (Solodov et al. 2007). We emphasize that the ignition energies mentioned 
above have been obtained under the strong assumptions of electron straight path, no 
beam divergence and without accounting for the self-generated EM fields. Atzeni et 
al. estimated that the scattering effects raise the ignition energy by about 20% (Atzeni 
et al. 2009). Thus, even for the ideal conditions assumed here, electron beam energies 
of several tens of kJ are needed to ignite a target. Similar results have been obtained 
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for the all-DT target design proposed for the HiPER project (Atzeni et al. 2008). 

5.2. Ignition energies of divergent electron beams with an assumed initial distribution 
function 

More realistic calculations can be performed by means of the partially integrated model 
mentioned above. In this model, the beam parameters are estimated from experiments 
or PIC simulations conducted at laser intensities and/or pulse durations lower than 
those required for FI. The main features of the relativistic electron source considered so 
far in FI simulations can be summarized as follows: 

i) Energy spectrum: It is normally asummed that the electron spectrum is 
given by the exponential distribution obtained in PIC simulations for sub-ps pulses. 
Relativistic Maxwellian spectra have been used also. Both distributions depend on the 
fast electron temperature or electron mean energy (E), typically fitted to experiments. 
For laser intensities <10 19 Wcm -2 , the so-called Beg's law (E) = 150 (In X 2 ) 1 ^ 3 is used, 
where (E) is in keV, In is the laser intensity in units of 10 17 Wcm -2 and A the laser 
wavelength in /im (Beg et al. 1997). For intensities around 10 19 Wcm~ 2 or higher, 
the ponderomotive scaling (Wilks et al. 1992) reproduces well experiments and PIC 
simulations. However, it gives electron energies well over the desired values around 2 
MeV for the laser intensities typical of the FI regime (>10 20 Wcm -2 ). 

Chrisman et al. (Chrisman et al. 2008), Haines et al. (Haines et al. 2009) and 
Kluge et al. (Kluge et al. 2011) have recently reported scaling laws that provide electron 
energies lower than those obtained by the ponderomotive scaling and of the same order 
than those predicted by the Beg's law. In principle, this is very favourable for FI 
because the optimal electron range is about 1.2 g/cm 2 (Atzeni 1999), which corresponds 
to an electron energy lower than 2 MeV (Solodov & Betti 2008) and a laser intensity 
about 2.4xl0 20 Wcm -2 assuming the Beg's scaling. Unfortunately, as it is discussed in 
Section 5.3, recent 3D PIC simulations have revealed that the mean energy of relativistic 
electrons is similar to that given by the ponderomotive scaling and their spectrum differs 
substantially from the exponential or relativistic Maxwellian distributions mentioned 
above (Strozzi et al. 2012). 

ii) Beam divergence: In most of the fast electron calculations carried out so 
far it has been assumed that the electron divergence is given by the beam effective 
propagation angle measured in the experiments. However, as electrons propagate in 
metals or plastics, whose resistivity is several orders of magnitude higher than that of 
the DT fuel, resistive collimation effects can be important (Bell & Kingham 2003). In 
this case, the initial fast electron divergence turns out to be substantially higher than the 
effective propagation angle measured (Stephens et al. 2004, Green et al. 2008, Norreys 
et al. 2009). For instance, to reproduce the full propagation angle of 35° found in the 
experiments conducted by Green et al. (Green et al. 2008), an initial electron divergence 
half-angle as large as 50° has to be assumed in hybrid calculations (Honrubia & Meyer- 
ter-Vehn 2009). Recent PIC simulations have shown initial divergence half-angles of 
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50-55° for FI conditions (Strozzi et al. 2012). Thus, guiding mechanisms should be 
envisioned for these highly divergent beams in order to have a good coupling efficiency 
with the dense fuel. 

It is also important to account for the dependence of the divergence angle on the 
electron kinetic energy observed in PIC simulations. It can be taken into account 
in a simplified fashion by assuming that the initial divergence angle is given by the 

1 /9 

ponderomotive scaling tan# = [2/(7 — 1)] ' (Quesnel & Mora 1998), where 9 is the 
polar angle and 7 the relativistic Lorentz factor. It is conjectured that an electron of 
energy (7 — l)m e c 2 is emitted with a divergence half-angle randomly selected between 
and 6. Thus, high energy electrons are well collimated while low energy electrons have 
an almost isotropic distribution. This dependence on electron energy is important to get 
a reasonable high energy coupling between the beam and the dense fuel. However, recent 
3D PIC simulations have shown that the energy spectrum and the angular distribution 
of the fast electron source are independent of each other, i.e. the divergence angle is the 
same for all electrons (Strozzi et al. 2012). This implies a strong increase of the ignition 
energies when compared with those obtained from the ponderomotive scaling. 

iii) Beam radius. Experiments and PIC simulations of relativistic LPI with foil 
targets show that the size of the fast electron beam is greater than that of the laser 
beam (Stephens et al. 2004). This effect has to be taken into account in fast electron 
transport calculations (Honrubia et al. 2006). In the case of electron acceleration in 
re-entrant cones, it is found that the radius of the beam is approximately equal to the 
cone outer radius. This is also true in the double cones described in Section 6.7, where 
the vacuum layer between the two cone walls force fast electron propagation towards 
the cone tip (Nakamura et al. 2007, Cai et al. 2010, Debayle, Honrubia, d'Humieres & 
Tikhonchuk 20106, Debayle, Honrubia, d'Humieres & Tikhonchuk 2010a). 

iv) Conversion efficiency. The laser-to-fast electron conversion efficiency 
obtained in PIC simulations ranges from 30 to 50% for the electron energies relevant 
for FI, e.g. E > 250 keV (Debayle, Honrubia, d'Humieres, Tikhonchuk, Micheau & 
Geissler 2010, Strozzi et al. 2012). 

As an example of the ignition calculations assuming an initial fast electron 
distribution carried out so far, let us discuss the ignition energies of the idealized DT fuel 
configuration shown in Fig. 5.2(a) (Honrubia & Meyer-ter-Vehn 2009). It is assumed 
that the fast electron energy and divergence are given by the ponderomotive formulas 
with A = 0.527 /zm multiplied by a scale factor. Despite it is technologically challenging, 
the 2 nd harmonic of the Nd laser has been considered to reduce the electron energy 
(Solodov et al. 2008, Honrubia & Meyer-ter-Vehn 2009, Strozzi et al. 2012). The DT 
fuel has an initial super-Gaussian density distribution, 498exp(— (i?/45) 4 ) gcm~ 3 where 
R is the distance to the centre in /iin, sited on a density pedestal of 2 gcm~ 3 . The total 
DT mass is 0.39 mg and the initial temperature 300 eV. As the cone tip is not included in 
the simulation box, its effect on fast electron transport (beam filamentation, scattering 
and energy loss) is accounted for indirectly via the initial fast electron distribution 
function. Calculations have been performed with the hybrid code PETRA for fast 
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electron transport (Honrubia et al. 2005) coupled to the radiation-hydrodynamics code 
SARA (Honrubia 1993) run in 2D Eulerian mode and cylindrical r-z geometry. 

Fast electron energy deposition takes place via Ohmic heating due to return currents 
and classical Coulomb scattering (collisional drag). Ohmic heating is important only 
in the low-density DT plasma, while collective behaviour is suppressed and energy 
deposition takes place almost exclusively by collisional drag in the dense core (Solodov 
et al. 2008, Honrubia & Meyer-ter-Vehn 2009, Johzaki, Nakao & Mima 2009, Strozzi 
et al. 2012). It is important to emphasize that collective effects may play a major 
role for core heating, but in an indirect way: self-generated B-fields may collimate the 
relativistic beam improving the coupling efficiency substantially. This effect appears 
to be very beneficial for FI: without resistive collimation there is little hope to ignite 
a pre-compressed target with reasonable beam energies due to the high divergence. 
However, as the beam collimation decreases strongly with the electron divergence angle, 
Eq.(9) (Bell & Kingham 2003), its importance in the FI scenario will depend on the 
full characterization of the fast electron source, which is not possible today neither 
by experiments nor PIC simulations. In addition to the beam collimation, resistive 
filamentation is the other collective effect that can play a role in electron-driven FI. It 
has been observed in the simulations of the ideal target shown in Fig. 16(a) (Honrubia 
& Meyer-ter-Vehn 2006), the imploded targets of Ref. (Solodov et al. 2008) and the 
experiments with solid and compressed plastic targets analyzed in (Solodov et al. 2009). 
It was shown also in Fokker-Planck simulations of sub-ignition targets (Johzaki, Nakao 
& Mima 2009). Its effects on fuel ignition have not been studied in detail yet. 




50 100 

z[(im] 



® 15- 

150 20 



divergence half-angle [°] 



Figure 16. Left: (a) Initial target density in gcm~ 3 . The halo surrounding the 
dense core has a density of 2 gcm~ 3 . Right : (b) Electron beam ignition energies as 
a function of the divergence half-angle and the electron mean energy (E) . The beam 
parameters for the case (35°, 1.6 MeV) are 20 fim radius (HWHM) and 18 ps pulse 
duration (FWHM) (Honrubia & Meyer-ter-Vehn 2009) 
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The minimum ignition energies of the target of Fig. 16(a) as a function of the 
initial divergence half-angle with the mean energy as a parameter are shown in Fig. 
16(b). For small divergences, the beam is strongly collimated and the ignition energies 
are even lower than those obtained in Section 5.1 for perfectly collimated beams. For 
higher divergences, the beam collimation is less important and ignition energies increase 
more than proportionally with the divergence half-angle. It is worth noting that, even 
for the ideal case shown in Fig. 16, to get ignition with a 50 kJ electron beam requires 
a substantial reduction of the divergence half-angle, from the 50-55° obtained from 
experiments and PIC simulations to the 35° shown in Fig. 16(b) for 2 MeV electrons. 

Similar calculations have been reported by Solodov et al. for the direct-drive capsule 
mentioned in Section 5.1 (Solodov et al. 2007). Simulations have been performed with 
the hybrid-PIC code LSP (Welch et al. 2001) coupled to the radiation- hydrodynamics 
code DRACO (Radha et al. 2005). A relativistic electron beam is injected 125 /xm 
from the target centre at a time when the maximum DT density is slightly above 500 
g/cm 3 . Assuming that beam electrons have a relativistic Maxwellian distribution with 
a mean energy of 2 MeV and a divergence half-angle of 20° (HWHM), Solodov et al. 
found a minimum ignition energy about 43 kJ (Solodov et al. 2008). They also found an 
important resistive collimation of the fast electron beam. The ignition energy increases 
strongly with the divergence half-angle 6, being 63 and 105 kJ for 9 = 30° and 40°, 
respectively. These energies are higher than those shown in Fig. 16(b) due to the higher 
stand off distance between the electron injection surface and the target centre (125 vs. 
100 /im) and the higher areal density of the coronal plasma surrounding the core of the 
imploded target (Solodov et al. 2008). For instance, Fig. 16(b) shows that the ignition 
energy for a beam with (E) = 2 MeV and 6= 35° is 50 kJ while it is around 80 kJ 
for the target of Ref. (Solodov et al. 2008). Assuming in this last case a laser-to-fast 
electron conversion efficiency of 40%, even with the strong assumptions of no cone tip 
and a divergence half-angle of 35° only, a multi-PW laser with more than 200 kJ would 
be needed to ignite such a target. 

The calculations presented in Refs. (Solodov et al. 2008, Honrubia & Meyer- 
ter-Vehn 2009) are quite ideal because fast electrons are injected just on the coronal 
plasma surrounding the fuel assembly. Indeed, fast electrons have to pass through the 
overcritical plasma sited inside the cone, if any, and always through the cone tip before 
reaching the dense fuel. Johzaki et al.(Johzaki, Nakao & Mima 2009, Johzaki et al. 2009) 
have studied the role played by the cone tip in the fast electron transport towards the 
dense core. They conclude that high-Z materials deteriorate substantially the quality 
of the electron beam even for cone tips as thin as 10 /xm. This is due to the collisional 
drag with plasma electrons and scattering with ions, resulting in an important reduction 
of the energy coupling of the fast electron beam to the dense core. The use of lighter 
materials in the cone tip, such as CH or DLC (Diamond-Like Carbon), mitigates drag 
and scattering effects. These light materials produce a manageable degradation of the 
beam and, at the same time, can support the shock and the jet coming towards the 
cone at the shell collapse time. DLC is preferred to the CH for its higher density. 
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The effect of the cone tip material in sub-ignition targets has been quantified by 
Johzaki et al. within the context of the FIREX-I project (Johzaki, Nakao & Mima 2009). 
The simplified target configuration considered in this project is a CD plasma with a 
Gaussian density profile in radius (peak density = 200 gem -3 , FWHM diameter = 20 
/an) placed on a density pedestal of 0.2 gem -3 . The target areal density is 0.2 gem -2 and 
the initial temperature is set to 400 eV. This imploded target configuration is heated 
during 5 ps by a 5 kJ, 1 MeV fast electron beam with a Gaussian radial profile of 30 /xm 
diameter (FWHM). The beam divergence full angle is 20° and the electrons are injected 
50 /im away from the core centre. As the electron beam size is comparable to the 
imploded core size and the electron range is higher than the target areal density, a hot 
spot is not generated and, instead, the whole core is heated. Johzaki et al. have shown a 
reduction of about 50% in the core heating energy and 60% in the peak ion temperature 
when a 10 /im thickness gold tip is present. This reduction is much lower, 10% and 
11%, respectively, for lighter materials such as CH (Johzaki, Nakao & Mima 2009). 

Another issue that can affect dramatically the fast electron beam coupling with the 
dense core is the existence of plasma inside the cone. As has been reported by Johzaki 
et al. (Johzaki et al. 2011), the ASE (amplified spontaneous emission) laser pre-pulse 
generates a long scale, low density plasma inside the cone. As a result, fast electrons are 
generated by the interaction of the main pulse with the pre-plasma relatively far from the 
cone tip and with much higher energy than without pre-plasma. This has been confirmed 
by PIC simulations (Johzaki et al. 2011), which have shown harder electron spectra 
and a substantial reduction of the number of electrons with energies lower than 5 MeV 
when a pre-plasma is present. Too hot electrons reduce substantially the electron energy 
deposition in the core and thus the coupling efficiency Integrated FI experiments carried 
out at the ILE have evidenced this effect, observing a coupling efficiency substantially 
higher when the laser pre-pulse energy is reduced (Shiraga et al. 2011) . 

An improvement of the standard field calculation model used in fast electron 
transport calculations has been proposed recently by Johzaki et al. (Johzaki, Nakao 
& Mima 2009) and Nicolai et al. (Nicolai et al. 2011). They have pointed out that the 
generalized Ohm's law and, in particular, the pressure gradient term, which yields an 
azimuthal B-field proportional to VT e x Vn e , can play a role in transport calculations. 
This term is important when the directions of the gradients of electron temperature, 
VT e , and electron density, Vn e , are not parallel, as occurs at the outer regions of the 
core (see Fig. 1(a)), where the electron temperature gradient is directed toward the 
beam axis and the density gradient is directed toward the core centre. The B-field 
due to VT e x Vn e has a direction opposite to the resistive collimating B-field mentioned 
above and its main effect is to scatter electrons away from the core. However, the growth 
of this B-field is relatively slow becoming important after several ps. For instance, the 
energy deposition in the target depicted in Fig. 16(a) heated by the beam defined by 
the parameters 9 = 35° and (E) = 1.6 MeV is reduced by 23% when the generalized 
Ohm's law is taken into account (Nicolai et al. 2011). A similar reduction (21%) has 
been found in the calculations for the simplified FIREX-I target discussed above when 
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the term VT e x Vn e is included (Johzaki, Nakao & Mima 2009). 

5.3. Ignition calculations with a PIC-based electron source 

A first characterization of the fast electron source in the FI scenario via 3D PIC 
simulations has been reported recently (Strozzi et al. 2011, Strozzi et al. 2012, Kemp & 
Divol 2012). One of the main conclusions of this study is that the initial distribution 
of fast electrons can be factorized as the product of two independent functions of angle 
and energy. The angular distribution is super-Gaussian, exp{— (9/A9) }, with A9 = 
90° and a mean divergence half- angle (8) = 52°. This high divergence can be explained 
by the curved geometry of the electron acceleration region (Schmitz et al. 2012) and 
by the electron scattering by the oscillating magnetic field generated by the Weibel 
instability close to the cut-off surface (Adam et al. 2006, Debayle, Honrubia, d'Humieres 
& Tikhonchuk 20106). The energy spectrum can be fitted by a quasi two-temperature 
profile. The first component is due to electron acceleration near the cut-off surface and 
has a temperature substantially lower than that given by the ponderomotive scaling T p , 
while the second component is due to electron acceleration in the subcritical plasma and 
has a temperature higher than T p . The overall electron mean energy is slightly higher 
than the ponderomotive temperature T p , but only 24% of the injected energy is carried 
by electrons with energies lower than T p . The main differences with the electron source 
assumed in Section 5.2 are the much higher electron energies and divergences, and the 
energy independent angular spectrum. 

Strozzi et al. (Strozzi et al. 2011, Strozzi et al. 2012) have performed integrated 
simulations assuming the PIC-based electron source presented above. The target used is 
shown in Fig. 2(a), where the DT fuel has initially a super-Gaussian density distribution, 
440exp{— (R/70) } gcm~ 3 where R is the distance to the centre in /im, sited on a 
background DT plasma of 10 gem -3 . A DLC cone is included in the simulation box. 
The DT mass is 0.57 mg and the initial temperature is set to 100 eV. An electron 
beam with quasi-uniform radial profile, exp{— In 2 (r/rb eam ) } with rb eam = 18 /zm, and 
constant intensity in time (from 0.5 to 18.5 ps, starting from a linear ramp from to 
0.5 ps) impinges on the target. The transport and energy deposition of fast electrons 
from the cone to the core has been simulated with the hybrid-PIC code Zuma (Larson 
et al. 2010) coupled to the radiation-hydrodynamic code Hydra (Marinak et al. 2001). 
The Zuma code includes the generalized Ohm's law. The details are given in (Strozzi 
et al. 2011, Strozzi et al. 2012). 
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Figure 17. (a) Initial target density in gem -3 . The red line indicates the source 
cylinder where fast electrons are injected. The 8 gem -3 (= 2.3 x solid diamond) 
carbon cone is coloured in blue for clarity, (b) Fusion yield vs. total injected electron 
energy, for Zuma- Hydra runs with an artificially collimated source A6>= 10°. ri, eam = 
10 (im for black squares with solid line, 14 /jm for red circles, 18 /xm for blue triangles 
and 23 /xm for green crosses. The blue triangle with Ef ast = 132 kJ is the lowest value 
that can be deemed to be ignited (Strozzi et al. 2012). 



The fusion yield as a function of the fast electron energy for an artificially collimated 
electron source (A6= 10°) and for several values of r^ eam is shown in Fig. 17(b). It is 
worth noting that, despite the artificial source collimation, the minimum ignition energy 
is quite high, 132 kJ and roughly independent of the beam radius. This can be explained 
by taking into account that for large radii, the increase in the volume to be heated is 
balanced by the lower energy and penetration of fast electrons due to the lower laser 
intensity. Assuming the laser-to-fast electron conversion efficiency of 0.52 obtained in 
the PIC simulations, the 132 kJ electron beam corresponds to a laser mean intensity 
of about 1.4 x 10 21 Wcm -2 , which gives an electron mean energy of 8.2 MeV and a 
range (6.8 gem -2 ) greater than the fuel pL (6 gem -2 ) and much higher than the optimal 
deposition range for FI (1.2 gem -2 ). This is one of the reasons why the ignition energy 
is so high. We can conclude that, even if a collimated electron beam could be generated, 
the ignition energy would be unacceptably high because fast electrons are too hot. When 
the fast electron divergence is included, simulations show ignition energies higher than 
1 MJ (Strozzi et al. 2012), making FI unfeasible, at least for the 'conventional' scheme 
discussed in this Section. The major limitations of the calculations presented here arise 
from the short time of the 3D PIC simulations (~0.36 ps), much lower than the 10 - 
20 ps pulse durations typical of FI, and from the lack of comparison with experimental 
data close to the true FI conditions. 
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6. Concepts for Controlling Transport 

6.1. Motivation 

From the preceding sections it is clear that large fast electron divergence angles can 
severely reduce the fast electron to hot spot coupling efficiency, thus raising the required 
ignitor pulse energy to levels at which FI is no longer a viable ICF scheme. In the first 
instance this is just a consequence of ballistic transport. Suppose that all other processes 
can be neglected, and that coupling is therefore dominated by ballistic transport with 
ideal stopping in the hot spot. When the stand-off distance, D, is much larger than 
the source spot radius, t\l, one would therefore expect the coupling efficiency ,%, to be 
roughly equal to the hot spot area divided by the fast electron beam cross-section at 
the stand-off distance, which is, 

2 

X = 2 J D 2 (l-cos# 1/2 )' (90) 

If 9i/2 is large then \ is limited to \ ~ r L/-^ 2 - For fhs/D = 0.2 this is 5%. 
Conversely, Eq. 90 indicates that for r^jD = 0.2, one needs #1/2 = 16.3°. Assuming 
that the laser absorption physics cannot be easily engineered to produce such a low 
divergence angle, there is a clear need to seek additional means of controlling the 
transport of the fast electrons. 

In the first instance this might be 'self-collimation' of the fast electron beam through 
the resistively generated magnetic field that should be produced around the fast electron 
beam (see §3.2). The ability of this field to produce significant focussing or pinching 
of the fast electron beam is at least doubtful in light of experimental and simulation 
results that have been obtained in the last few years. 

If self-collimation cannot be relied upon, then the FI scheme has to be adapted in 
some way so as to ensure effective transport of the fast electrons by some other means 
(if at all possible). This might be possible through optical engineering of the laser pulse 
and a thorough understanding of the laser absorption process. Alternatively it might 
be possible to exploit fast electron transport physics, and it is this that this Section is 
concerned with. 

6.2. Self- Collimation of the Fast Electron Beam 

In §3 of this review, the criteria for a fast electron beam to 'self-collimate' as derived 
and studied in the work of Bell and Kingham. By this, we mean the collimation of a 
fast electron beam even in a a homogeneous background plasma due to the resistive 
magnetic field that is generated due to the curl of the fast electron current density, and 
thus the curl of the electric field. A simplified version of the Bell-Kingham criteria can 
be derived as follows. One starts with an estimate of the magnetic field that is resistively 
generated, 

d — « ™l (91) 

dt R ' l } 



B 



,5/21 2/5 



at + T° /z , (92) 



?i/ 2 tan 0i/2 



at + T 5/2 l 2/5 , (94) 
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where i? is the beam radiusj/ the fast electron current density,?] the resistivity, and £ 
is the time. If we use the Spitzer resistivity (77 = aT~ 3 ^ 2 ), then this can be integrated 
to obtain, 

3n e 

where a = 2ajj/3en e , and To is the initial temperature. Next one estimates the angle 
(0) that a fast electron will be deflected through over the distance that the unperturbed 
beam takes to double its radius, 

e= eB \ . (93) 

7m e c tan 6*1/2 

where 0i/2 is the divergence half-angle of the beam. If we use the fundamental Bell- 
Kingham criterion for collimation — that collimation occurs when = 0i/2 — then we 
can combine Eq.91 and 93 (and substitute power balance for jf) to obtain, 

3ene 
jm e c(3I L 

where /3 is the laser to fast electron conversion efficiency, II is the laser intensity, and 
e is the average fast electron energy. Inserting typical numbers into Eq.94 leads to the 
conclusion that self-collimation is only likely to happen for 0!/ 2 <20-30°. 

Calculations of fast electron transport relevant to ignition scale FI indicate that 
this is a reasonable estimate for the regime in which self-collimation is sufficient 
(e.g. (Honrubia & Meyer-ter-Vehn 2009) ;see Section 5). There is good evidence, 
however, that the actual fast electron divergence half-angle under ignition scale 
conditions will be significantly greater than 30°. Some of this evidence is experimental 
(Lancaster et al. 2007, Green et al. 2008) coupled with theoretical and numerical 
interpretation (Honrubia & Meyer-ter-Vehn 2009). Other evidence comes from large 
scale numerical simulations of laser absorption ,e.g. (Debayle, Honrubia, d'Humieres 
& Tikhonchuk 20106). With divergence half-angles that are slightly in excess of 30°, 
magnetic field generation is still highly beneficial in terms of improving the coupling. 
However the evidence suggests that the divergence half-angle could be in excess of 50°, 
and under these conditions self-collimation does little to prevent very poor coupling to 
the hot spot (for typical stand-off distances). 

6. 3. Resistive Guiding of Fast Electrons 

From the induction equation in the hybrid approximation, 

rfVK 

_ = ?7 Vxj / + Vr7XJ / , (95) 

one can see that just as there is a term which indicates that magnetic field grows to 
drive fast electrons into regions of higher fast electron current density (first term) , there 
is also a term that indicates that magnetic field grows to drive fast electrons into regions 
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of higher resistivity (second term) §. The 77 V x jf term is the effective responsible for 
the 'self-collimation' described in the preceding section. 

Resistive guiding exploits the second term, the Vi] x jf term (Robinson & 
Sherlock 2007). At sufficiently high temperatures, all materials will follow a Spitzer-like 
resistivity in which rj oc Z. Therefore if one structures a target by using two materials 
with different Z, in principle the higher Z material should confine and guide the fast 
electrons as magnetic fields are generated at the interace between the two materials 
where the V?7 x jf term will be large. 

It is relatively straightforward to see that the collimating fields generated by 
resistivity gradients must be at least as strong as those generated by the r/V x j/ term. 
If we denote the 77V x jj term by Bi, and Vi] x j f by B 2 , then we can see that the 
magnitudes are approximately, 

Bi « ^, (96) 

Kb 

B 2 » ^L, (97) 

where Rj, is the radius of the fast electron beam, and L int is the scale-length associated 
with the transition in resistivity. As it is possible to envisage transitions in resisitivity 
with scale-lengths, L int <C i?&, due to the sharp interfaces produced by target 
engineering, there will be range of circumstances in which resistivity gradients are 
capable of producing powerful confining magnetic fields. 

6.4- Preliminary Studies of Resistive Guiding 

In the work of Robinson and Sherlock (Robinson & Sherlock 2007) this concept was 
explored using a 'hybrid- VFP' code and very good guiding was demonstrated. This 
investigation concentrated on the conditions close to those of laboratory experiments 
involved solid density foils and PW class lasers that deliver a few hundred Joules in 
about a picosecond. The possibility that 'cold target' effects might cause a problem 
for the concept was also investigated and it was shown that a small temperature range 
over which the resisitivities are inverted can be tolerated. One might expect that once 
collimation has been occurred, that the 77V x jj will then act to reinforce this collimation 
and sustain it, even if the Vi] x jf greatly diminishes. Robinson and Sherlock showed 
some evidence for this sort of positive feedback in this concept. These initial simulations 
were done in 2D Cartesian coordinates, and only considered a guiding structure that 
was perfectly aligned with the axis along which the fast electrons were injected. In 
Fig. 18 results from one of these early simulations are shown, including the material 
composition of the target and the fast electron density after several hundred fs. 

§ This is Davies' qualitative interpretation of this equation (Davies 2002). 
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Figure 18. (Left) Plot of target Z, and (Right) Plot of fast electron density (login) 
after several hundred fs in an early simulation carried out using the 2D hybrid- VFP 
code, LEDA. 



Some promising results were obtained in early experiments carried out using PW 
class lasers and 'cold' solid-density targets (Kar et al. 2009, Ramakrishna et al. 2010), 
and more simulation work was carried out using a 3D particle hybrid code to help 
interpret these experiments. This combination of numerical studies and experiments 
gives one good reason to believe that resistive guiding is a genuinely realizable effect 
that might be exploited in FI. Subsequent work in this area has therefore focussed on 
two major concerns: (i) Can resistive guiding work under the conditions that FI imposes 
on it?, and (ii) How precisely can we exploit resistive guiding in a realistic FI scenario? 



6.5. Resistive Guiding at High Energy Scales 

Applicability of the fast-electron collimation scheme exploiting resistivity gradients 
to fast ignition depends on two factors: (i) The collimating magnetic fields need to 
persist during the entire ignition pulse; (ii) The high resistivity path must survive the 
compression resulting from the implosion. In an ignition pulse, the material along the 
path of electron propagation is heated to keV temperatures. At such high temperatures, 
the resistivity of such a material can become less than the resistivity of the surrounding 
plasma and the resistivity gradients can be inverted. The question is whether or not 
the inverted resistivity gradients can cause a magnetic field reversal from collimating to 
de-collimating. It can also be difficult to maintain the guiding structure extending to 
the dense fuel up to the time of significant compression of the fast-ignition target. This 
problem should be addressed either by developing target designs for which the damaging 
effect of the implosion is minimized or placing the guiding structure inside a protective 
cone as was suggested in (Robinson et al. 2012). 

Collimation of high-energy electron beams in the wire-like structures has been 
studied (Solodov, Betti, Anderson, Myatt, Theobald & Stoeckl 2010). Simulations 
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using the hybrid-PIC code LSP were performed for a 40 /zm diameter copper wire 
embedded in aluminum. A 10 ps (constant in time), 2 MeV mean-energy, relativistic- 
Maxwellian electron beam with divergence half- angle of 67°, and total energy of 20 kJ 
was injected into the wire. The beam was found to be effectively collimated for the whole 
duration of the electron pulse. About 65% of the injected electrons were collimated on 
the length of the collimating structure of 150 /xm. The resistivity gradient at the wire 
boundary was found to be inverted because of the wire heating by fast electrons in less 
than 0.5 ps after the beginning of the electron pulse. The collimating magnetic field, 
however, persisted because the magnetic field had two components: one generated by 
the resistivity gradients and the other by the return current gradients. Initially, the 
collimating magnetic field was generated by the resistivity gradients. The resulting 
collimation caused large current density gradients that enhanced the collimating field. 
The current density gradient offset the effect of the reversal of the resistivity gradient 
thus supporting a large saturated collimating magnetic field. Similar conclusions were 
obtained for lower-energy electron beams in (Ramakrishna et al. 2010). 

LSP simulations (Solodov, Betti, Anderson, Myatt, Theobald & Stoeckl 2010, 
Solodov et al. 2010) showed that high-energy electron beams can be guided by a mid- 
Z wire through the cone tip and coronal plasma of a fast-ignition target, subject to 
survivability of the wire during the implosion. The simulations utilized idealized cone- 
fuel configurations with and without a wire (Fig. 19). A 75 /im long, 40 /im diameter 
copper wire goes through the 25 /im thick tip of aluminum cone towards the pre- 
compressed deuterium fuel core. The core has a super-Gaussian density distribution 
400exp(— (r/50) 4 ) gcm~ 3 where r is the distance from the center in /mi, sited on a 
background deuterium plasma of 10 gem -3 . The initial temperature of 100 eV was 
assumed, ionization and radiative cooling were modeled for copper and aluminum. The 
simulations used a 40 kJ, 10 ps, 1.6 MeV-mean-energy relativistic Maxwellian electron 
beam with initial divergence half-angle of 55°(oc exp(— (6/6o) 2 ), with 6q =67°), constant 
temporal profile, and a super-Gaussian radial profile 400exp(— (r_|_/20) 4 ), where r± is 
the distance from the beam axis in /im, injected at the inner side of the cone tip. 
Comparison simulations were performed in which the cone and the wire were replaced 
by a single copper cone. Fast-electron energy deposition in the cylindrical region (see 
Fig. 19) with a diameter of 60 /im and a length of 40 /im (so-called ignition region) 
was calculated. The simulations confirm that the fast-electron coupling to the core is 
significantly improved with the wire: 45% coupling efficiency to the ignition region in 
the cone-wire case versus 7% without a wire. It can be seen in Fig. 19 that fast electrons 
are effectively collimated and guided by the self-generated resistive magnetic fields at the 
interface of the copper wire and surrounding lower-Z plasma. Collimation of electrons 
to the dense fuel in a wire-like structure has been also confirmed by hybrid simulations 
of J. Honrubia and D. Larson using codes PETRA and ZUMA (Solodov et al. 2010). 

The question of wire survivability was not addressed but it was noted that it may be 
difficult to maintain a clean high-resistivity path to the dense core at the time of peak 
compression in an imploded capsule. Detailed target design studies using radiation- 
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hydrodynamics codes are required to show if such a divergence-mitigating structure can 
be assembled in an actual implosion. 
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Figure 19. LSP simulations (Solodov, Betti, Anderson, Myatt, Theobald & 
Stoeckl 2010) predicting significant increase in energy coupling to the compressed core 
by the self-generated magnetic field at the copper wire interface. 



6.6. Advanced Uses of Resistive Guiding 

The use of resistive guiding in FI may not necessarily be restricted to the cone-wire 
scheme discussed above. Other schemes have been suggested that do not require 
placing an element outside of the cone. Instead, such schemes suggest putting guiding 
structures in an insert in the cone tip. This may have the advantage of being more 
robust with respect to the hydrodynamics of fuel assembly, although detailed radiation- 
hydrodynamics studies are required to confirm this. 

One, suggested by Schmitz (Schmitz et al. 2012), is to use a curved axisymmetric 
interface (ellipsoidal or paraboloidal) to produce an azimuthal magnetic field structure 
that acts as a curved mirror. A divergence fast electron spray will then have its angular 
spread reduced by the approximately specular deflection in the strong magnetic field 
which is localized at this interface, in a way that is analogous to a parabolic or elliptical 
mirror in ray optics. This comes at the expense (as in the case of an optical parabolic 
mirror) of increasing the radial extent of the beam. As the hot spot radius and laser 
spot radius are likely to be comparable in size, there will be limits on how much the 
radial extent of the beam can be increased to reduce angular spread (the diameter of 
the cone apex may also have to be limited for fuel assembly). This was not assessed in 
Schmitz's original proposal due to limitations on the size of the simulation domain. Some 
preliminary simulation results presented by Robinson suggest that some substantial 
benefits might still be obtained from the elliptical mirror approach (Robinson 2012) 
where a simple elliptical configuration yielded about an improvement in coupling into a 
target hot-spot region of about 2-3 (over an unguided case). 
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Another suggestion was made by Robinson and co-workers (Robinson et al. 2012) 
who suggested a 'magnetic switchyard' configuration: a series of concentric quasi- 
cylindrical guide shells immersed in a less resistive substrate. Strong azimuthal fields 
grow at the interfaces, confining fast electrons into the guide shells. Each guide shell 
will receive a portion of the fast electron population with limited angular spread about 
some mean angle. The guide shells then curve around in an arc which re-directs this 
mean angle to some distant region. As with the aforementioned mirror concept, this 
means that the switchyard must increase its radial extent beyond that of the source 
radius. The numerical simulations presented by Robinson showed that, at least for 
some configurations, an improvement in coupling of about 2-3 (also compared to an 
unguided case). 

To the same goal, Debayle et al. (Debayle et al. 2013) recently proposed a structured 
target made of narrow high- and low-Z elements of density decreasing in the axial 
direction. The magnetic modulations developing at the filament interfaces then decay 
away from the surface, leading to non-specular reflections of the fast electrons trapped 
inside the high-Z filaments. As a result, their local angular dispersion steadily decreases 
along their path. The capability of these targets to both guide the fast electrons and 
reduce their angular dispersion is, however, obtained at the cost of heavy constraints on 
target manufacturing. 

Finally, note that transverse resistivity gradients associated to density or 
temperature non-uniformities in the corona also have the potential to beneficially affect 
the fast electron transport. This was demonstrated both experimentally and numerically 
in Ref. (Perez et al. 2011) in the case of cylindrically-compressed foam targets. Insofar 
as they are injected shortly before the shock convergence, the fast electrons can be 
efficiently guided by the magnetic field growing at the shock front. On the other hand, 
non- uniformities can help drive resistive filamentation (Robinson et al. 2008), as well as 
the de-collimating effects noted in Section 5, so the exploitation of hydrodynamically 
induced non-uniformities requires careful examination. 

6. 7. Double Cone Approach 

In this concept, a vacuum gap is introduced in order to prevent those fast electrons 

travelling at large angles from escaping. The vacuum gap is introduced by employing 

a target using two concentric cones rather than one, hence 'Double Cone'. Due to the 

vacuum gap, the cone wall is isolated from the coronal plasma and the fast electrons are 

confined and guided to the tip by electrostatic and quasi-static magnetic fields formed 

in the vacuum gap region. 

The fast electron guiding using vacuum gap has been first proposed by Campbell, 
et al. (Campbell et al. 2003) , where a collimation of high energy electrons using planner 
plag/gap/foil structure was numerically shown and an idea to control of fast electron 
using a conical plag/gap/foil structure (see (Campbell et al. 2003)) was proposed for 
fast ignition application. This scheme can be applied for the beam guiding from cone 
tip to the core. 
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Figure 20. The double cone target employed in simulations by Cai (Cai et al. 2010). 



Nakamura et al. (Nakamura et al. 2007) have suggested a double cone target 
(Fig.20), where the vacuum gap is introduced into the side wall of cone target. On 
the basis of two-dimensional (2D) Particle-in-cell (PIC) simulation, they showed that 
the double-cone confines the electrons for hundreds of femtoseconds (fs) by the sheath 
electric field generated inside the vacuum-gap. However, the simulation time was limited 
to a few hundreds fs. Contrary to this, the core heating duration in practice is lOps 
order. So the reduction of sheath electric field inside the gap due to the plasma expansion 
from the cone wall in the early phase of core heating and then the failure of confinement 
was feared. Later, Cai et al. (Cai et al. 2010), carried out ps order 2D PIC simulations 
and demonstrated that the double cone is still effective in confining the high-energy 
electrons even for the gap width of a few microns (Fig. 21). 
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Figure 21. Natural Logarithm of electron energy density for (a) single cone, and (b) 
double cone at 1 ps in simulations performed in (Cai et al. 2010) 



After reduction of sheath electric field due to plasma expansion, the quasi-static 
magnetic field works to confine the fast electrons (Fig. 22). The quasi-static magnetic 
field has been generated due to a localized supply of high-energy electrons, originally 
produced at the inner-cone and the cone tip. This electron current coming from the cone 
tip produces a positive current inside the gap, while an opposite-directional surface- 
current is generated along the inner-surface of the outer-cone. The collaboration of 
these two currents generate a large quasi-static magnetic field inside the gap. These 
quasi-static fields continue to confine the high-energy electrons for longer than a few 
picoseconds. They showed that the double cone can reduce the beam energy loss from 
the side wall down to 1/3 of that for the single cone case. 
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Figure 22. PIC Simulations results from (Cai et al. 2010) showing : (a,b) time- 
averaged sheath electric fields at 330 fs and 1500 fs, and (c,d) time-averaged magnetic 
fields at 330 fs and 1500 fs. Fields are in normalized units, i.e. m e LJLc/e. 



Even if the double cone is used, the fast electron beam after penetrating the cone 
tip spatially diverges during propagation to the core due to its large divergence angle 
and then the enhancement of core heating rate may not be expected so much. In order 
to guide the fast electron beam close to the core, Johzaki, et al. (Johzaki et al. 2011) 
have proposed to extend the cone tip and vacuum gap (Fig. 23). They called it the 
extended double cone, which is a combination of the conical plag/gap/foil structure 
(Campbell et al. 2003) and the double cone (Nakamura et al. 2007, Cai et al. 2010). In 
this case, the fast electrons travel a long distance in the extended tip region, so that 
low-Z, but relatively dense material (e.g. DLC) was proposed as the tip material to 
reduce the collisional effects (Johzaki, Nakao & Mima 2009, Johzaki et al. 2009) . It 
was shown from the 2D PIC and FP simulations for core heating that the extended 
double cone with a short inner wall enhances the core heating rate more than four times 
when compared with the single wall cone case (see fig. 23 and table 1). 
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Figure 23. Schematic view of the extended double cone proposed in (Johzaki 
et al. 2011). 



The preliminary experiments have been conducted by Sakawa et al. (Sakawa 
et al. 2010) to prove the vacuum-shielding effect using an Al-Cu double-foil targets with 
and without a vacuum gap. The enhancement of the number of electrons detected in the 
target surface direction has been observed for the case with the vacuum gap compared 
to the case without vacuum gap, which demonstrates the fast electron confinement by 
the vacuum gap. 

For more realistic evaluations of the possibilities for the extended double cone, 
hydrodynamic modeling of the cones ignition-time structural distortion must be included 
in the simulations. Also, the integrated experiments are indispensable to prove the 
guiding performance. 



6.8. Axial Magnetic Field Approach 

There is another approach to divergence mitigation which employs an imposed magnetic 
field, rather than a field self-generated by the fast electron current. There are two main 
methods which have been suggested for imposing the needed multi-MG fields: flux 
compression in the fuel assembly implosion (Shay et al. 2012, Tabak et al. 2010), and 
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Cone Type 


Pre-Plasma 


^fe(%) 


J ?t^fe<ioMeV( % ) 


J ?fe^core( % ) 


r /L^core( % ) 


(Ti) DD (keV) 


Single Cone 


None 


18 (48) 


15 (39) 


16 


7.5 


0.75 


Single Cone 


\ p =10/im 


14 (36) 


4(11) 


4.7 


1.7 


0.35 


Ex. Dbl. Cone 


None 


31 


27 


62 


19 


1.27 


Ex. Dbl. Cone 


\ p =10/im 


20 


11 


28 


5.5 


0.7 


E.D.C. with short inner wall 


none 


41 


35 


79 


32 


1.58 


E.D.C. with short inner wall 


\ p =10/im 


23 


14 


36 


8.1 


1.01 



Table 1. Summary of heating results from Extended Double Cone calculations 
(Johzakiet al. 2011). 



laser-driven coils (Daido et al. 1986, Fujioka et al. 2012, Fujioka et al. 2013). A body of 

simulation work has been carried out at LLNL on the assembly of such fields, and the 

characterization of their advantages for electron transport, which we shall review here. 

The purpose of an imposed field is to spatially confine the fast electrons to small 

radius (perpendicular to the axial direction), and enhance their flux on the fuel. We 

distinguish confinement, or limiting the fast electrons from spreading in space, from 

collimation, or reducing their velocity-space divergence. A confining magnetic field 

generally will not collimate, so fast electrons emerge from a confining magnetic field with 

their original divergence. An estimate of the product of field strength times path length 

needed to confine a fast electron of velocity v is given by (Robinson & Sherlock 2007) 

m c 
BL > Kjptl - cos 9) K=^ = 17.0 MG • fim. (98) 

e 

(3 = |v|/c, 7 = (1 — /3 2 )~ 1//2 , and 9 is the angle between v and the z axis. For instance, 

a 3 MeV electron with 9 = 45° requires BL > 33.9 MG-/im to be confined. Keeping the 

field thickness smaller than the source spot size imposes L < 10 /an, or B > 3 MG. 

Flux compression (Velikovich 2012) exploits the frozen-in law of MHD, which states 
that the magnetic flux oc B-rfa enclosed by a good conductor of area a remains constant. 
As a decreases the field strength rises. A "good" conductor is one for which the resistive 
diffusion time ~ fi aL 2 of the magnetic field is much longer than the implosion time 
(a is the conductivity and L a field length scale). Implosions at the Omega laser have 
compressed axial seed fields of <0.1 MG to 20-40 MG in cylindrical (Knauer et al. 2010) 
and spherical (Chang et al. 2011, Hohenberger et al. 2012) geometry. 

Preliminary fast ignition implosion simulations with an initial seed magnetic field 
have been performed, using the MHD capabilities of Hydra (Shay et al. 2012) and 
Lasnex (Tabak et al. 2010). These have considered radiation-driven implosions around 
a re-entrant cone. Figure 24 shows the magnetic field in a Hydra simulation done by 
H. Shay, starting with a uniform, axial field of 0.1 MG. It is similar to the implosions 
presented in Ref. (Shay et al. 2012), and entails a radiation-driven beryllium ablator, a 
DT ice layer, and a carbon-tipped gold cone. The plot is taken at the time of peak fuel 
compression, and shows a field of 500 MG in the compressed fuel. However, the field 
does not diffuse far into the cone or its interior, where the field in the critical-density 
plasma is < 20 MG. The fast electrons will therefore be generated by short-pulse LPI 
in a region of relatively low field. This poses two challenges. First, electrons which 
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Figure 24. Magnetic field from a Hydra simulation with an initial axial field of 0.1 
MG, performed by H. Shay and detailed in the text. The solid or dashed black contours 
are material interfaces, the colored contours are density (the red, bark blue, and yellow 
values are indicated), and the solid background color is |B| on a log scale. Text labels 
indicate material. 



encounter a field that increases along field lines (e.g., an axial field that increases in the 
axial direction) are subject to magnetic mirroring and reflection. In addition, there is a 
standoff distance before the fast electrons reach a field strong enough to confine them. 

Imposed multi-MG axial fields have been studied using the coupled hybrid-PIC 
code Zuma and the rad-hydro code Hydra in Ref. (Strozzi et al. 2012) (see Section 5, 
for previous description). This work considered an idealized, spherical DT fuel assembly 
with a carbon cone, and simple initial field profiles. Fusion yields are presented in Fig. 
25. The electron source had a substantial divergence based on full-PIC LPI simulations: 
dN/ain oc exp[-(6»/A6») 4 ] with A0 = 90°, giving (9) = 52°; Q is the velocity-space solid 
angle element. A 50 MG uniform, initial axial field performed slightly better than an 
artificially collimated source with A9 = 10° (or (8) = 7°) with no imposed field. Both 
cases required ~ 130 kJ of fast electrons to ignite. A field of 30 MG needed almost twice 
as much energy to ignite, while a 10 MG field performed significantly worse (although 
still much better than with no initial field). This work explored the degradation due to 
mirroring in non-uniform field profiles, and presented the hollow magnetic pipe as one 
way to avoid mirroring. The effects on yield are presented in Fig. 25. The non-uniform 
field cases labeled BZ30-75, BZ50-75, and BZO-50 demonstrate the reduced coupling 
due to mirroring. The magnetic pipe, case BZ50-pipe, couples almost as well as the 
uniform field case BZ50. Figure 26 shows the pipe field envelope; the thick pipe (solid 
black outline) is used in Fig. 25. More work on assembling pipe field structures in 
fast-ignition implosions is needed to make this scheme viable. 

The role of different terms in Ohm's law has been studied in several recent works, 
such as (Johzaki, Nakao & Mima 2009, Nicolai et al. 2011, Strozzi et al. 2011, Strozzi 
et al. 2012) (models based on Ohm's law are detailed in Section 4.2 of the present 
review). All of them have shown that the fast electron coupling to a spherical fuel 
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Figure 25. Reprinted from (Strozzi et al. 2012): fusion yield vs. added fast electron 
energy for Zuma-Hydra runs with PIC-based divergence A9 = 90°, except for one with 
an artificially collimated source A8 = 10° and no initial B field (labeled DQ10). Other 
cases had no initial B z (DQ90), initial uniform B z = 10, 30, 50 MG (BZ10, BZ30, 
BZ50), non-uniform initial B z rising gradually from 30 to 75 MG (BZ30-75), gradually 
from 50 to 75 MG (BZ50-75), rapidly from 0.1 to 50 MG (BZ0-50), and the magnetic 
pipe with 50 MG peak field (BZ50-pipe). 



region is degraded with an extensive Ohm's law that includes terms beyond the resistive 
i]J e , especially Vp e . This produces an azimuthal Vn c x VT e magnetic field at the fuel 
interface that pushes fast electrons to large radius and away from the hot spot. 

Imposed magnetic pipes of differing orientation illustrate the potential benefits 
of axial (B z ) and azimuthal (.£?</,) fields. Both the field direction (z vs. <fi) and sign 
significantly affect its confinement properties, as do non-resistive terms in Ohm's law. 
Electrons are confined in radius by the radial component of the v x B Lorentz force, 
which is independent of v r and B r . For a simplified discussion we imagine r and <fi to 
be Cartesian, and neglect the change in v r and v^ due to free motion. This is valid for 
sufficiently small Larmor radius. The expected confinement based on electron orbits is 
as follows. Each sign of B z should confine electrons with one sign of v $ to a small radial 
excursion, and the other sign of v^ to a larger one. Since the electron source is expected 
to be uniform in v^, we expect comparable confinement from either sign of B z . For 
Bq < 0, however, both signs of v^ are well confined, while both are poorly confined for 
B$ > 0. We therefore expect B^ < to confine the best, B^ > to confine the worst, 
and both signs of B z to be intermediate and comparable to each other. 

Zuma-Hydra simulations with the same profile of initial |B| have been performed 
for four cases all with peak magnitudes of 50 MG: B z > 0, B z < 0, B^ > 0, and 
B^ < 0. For the B z cases the magnetic field is found from a vector potential A^, so 
that B r is included to satisfy V ■ B = 0; no such B z or B r is needed for B^,(r,z). 
These calculations show a different ordering of confinement quality than the simple 
orbit discussion. Figure 26 plots the fusion yield for various initial fields. As expected 
from orbits, Ba, < performs the best, and B^ > the worst. Although both B z have 
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Figure 26. Left: Initial magnetic field envelope \B Z \ or \B^\ for integrated Zuma- 
Hydra runs with thick and thin "pipe" . Right: fusion yield vs. injected fast electron 
energy for thick pipe with B z q > (taken from Ref. (Strozzi et al. 2012)) and thin 
pipes with different B signs and directions. 
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Figure 27. Fast electron power reaching right edge in simplified Zuma-only 
simulations with initial, "thin" magnetic pipes as in Fig. 26. 



intermediate performance, B z < confines better than B z > 0. The results presented 
in Ref. (Strozzi et al. 2012) unfortunately used B z > 0, and would be better with the 
opposite choice. Moreover, the better confinement with an imposed B^ is promising for 
self-generated field approaches, which usually give rise to a 5^. 

Simplified Zuma-only runs (no coupling to Hydra) were performed to study 
magnetic pipes of different field orientation. They indicate that the different coupling 
for the two signs of B z results from different magnetic field evolution, which occurs only 
when non-resistive terms are included in Ohm's law. These runs used a uniform DT 
plasma of p = 10 g/cm 3 and 100 eV, and a divergent electron source with A8 = 90°. 
Figure 27 depicts the power of fast electrons reaching the right edge, inside the pipe 
radius, and counting only 1.3 MeV of kinetic energy per electron (the maximum 
deposited in a DT hot spot of optimal depth pAz = 1.2 g/cm 3 ). For both the left 
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panel (E = 0, B = constant) and center panel (E = i]3 e , B evolves), the ordering is as 
expected from orbits. However, when a more extensive Ohm's law is used, the two signs 
of B z perform differently, with B z < coupling better than B z > 0. Work is ongoing to 
elucidate the difference in field dynamics. 

Magnetic field generation by laser-driven coils (Daido et al. 1986) was demonstrated 
experimentally in the 1980s, and has recently been suggested as a divergence mitigation 
approach for fast ignition (Fujioka et al. 2012). In this scheme, a coil connects two 
plates, one of which is struck by a kJ-class, long pulse (~ ns) laser. This generates hot 
(~ 10 keV) electrons by resonance absorption, which reach the other plate and set up a 
potential difference between the two plates. A large, transient current flows through 
the coil and induces a large magnetic field. Experiments at GEKKO-XII (Fujioka 
et al. 2012) have recently made fields of 10 MG. Ref. (Courtois et al. 2005) reports 
smaller-scale experiments at the VULCAN laser which generated fields of 0.1 MG, as 
well as a simple model of the system. This approach eliminates the need to create 
large fields in an implosion, and allows one to consider novel field configurations. For 
instance, the field could peak on the short-pulse laser side of the re-entrant cone, so that 
the mirroring effect pushes fast electrons toward the fuel. The field is generated over 
a few ns, which is only a fraction of the duration of a typical fast ignition implosion. 
It can therefore be timed such that the implosion modifies the field in a limited way. 
Realistic rad-hydro-MHD modeling is needed to validate the laser-driven coil approach 
in specific geometries, and when coupled to specific implosions. 

7. Conclusions 

Although Fast Ignition ICF was proposed nearly 18 years ago, in-depth simulation 
studies of the fast electron transport aspect of the problem have only really been done 
over the past 7 years or so (Section 5). In this review we have looked at the following 
aspects : 

(i) Basic Physics: The fundamental physics of fast electron transport is thought 
to be well understood. Where doubts do exist is in the ability to calculate the properties 
of dense matter, and collisional processes to very high accuracy. Another aspect 
which hasn't been fully explored is the extent to which kinetic micro-instabilties affect 
transport, and whether significant corrections could be included in hybrid models. These 
aspects need to be addressed, but it is thought that these are unlikely to radically alter 
the outlook for Fast Ignition. 

(ii) Simulation Methods: In Section 4 it was shown that there are now a number 
of different simulations models and a large number of codes. These are nearly all 
extendable models that can include a wide variety of physics, and most have some 
"hybrid" character. There are a number of codes that include hydrodynamic motion, 
radiation transport, and thermonuclear burn - - representing a true merging of the 
hybrid model with standard ICF rad-hydro models. Thus there are codes that can 
tackle most of the full problem. This is subject to the inclusion of a realistic fast electron 
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source and the approximations inherent in "hybridization" . Nonetheless, as this review 
shows, this is sufficient to use computer simulation to make a first assessment of the 
viability of the Fast Ignition concept. 

(iii) Ignition-scale Calculations : In recent years there have been a number 
of studies (see Section 5) which have attempted to use the aforementioned numerical 
models to assess the viability of Fast Ignition at (or close to) ignition scale. This has 
ranged from rather idealized calculations through to calculations which have attempted 
a good degree of realism. These calculations have shown that achieving modest or 
"attractive" ignition energies (<100 kJ of laser energy) is difficult. The good laser-to- 
fast electron conversion efficiencies that are seen in certain PIC calculations (40-50%) 
are certainly beneficial, but the degradation of the coupling into the hot spot due to 
high fast electron divergence is considerable. 

(iv) Controlling Transport : As the divergence problem is a serious one, a 
number of studies have looked at modifications to the FI scheme that will allow the 
flow of fast electrons to be sufficiently controlled as to improve the coupling into the 
hot spot. In Section 6 a number of interesting possibilities were discussed including 
the exploitation of resistivity gradients, axial magnetic fields, and double cones. This 
is a highly active area, which may eventually produce an attractive solution to the 
divergence problem. 

The problem of simulating fast electron transport in Fast Igntion is not a trivial 
problem. Nontheless, as this review shows, there has been considerable progress over 
the last seven years or so. The main challenge that is now being faced is how one can 
modify transport so that relatively modest ignition energies can be achieved given a 
divergent fast electron source. A number of concepts have been proposed, and thus it 
is still entirely possible that an attractive Fast Igntion point design will emerge in the 
coming years. 
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